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Abstract 


In  a  review  of  studies  on  the  seismic  phase  Lg,  we  describe  its 
particle  motion,  dispersion,  spectral  content,  mode  of  propagation, 
and  magnitude-scale;  we  also  tabulate  the  regional  velocity,  atten¬ 
uation,  and  propagation  efficiency  for  this  seismic  phase. 

The  characteristics  of  Lg-wave  propagation  in  the  eastern  and 
western  United  States  are  compared  with  those  in  different  regions 
of  the  Soviet  Union.  Possible  discriminants  such  as  (i)  Lg  vs.  P 
amplitudes,  (ii)  Lg/P  amplitude  ratios  as  a  function  of  distance, 

(iii)  group  velocities  of  Lg  at  amplitude  maxima,  and  (iv)  Lg  energy 
ratios  are  found,  similar  to  attenuation  and  group  velocity,  to  be 
highly  deperfdent  on  the  propagation  path.  The  valid  application  of 
these  quantities  to  the  problem  of  earthquake-explosion  discrimina¬ 
tion  will  therefore  require  regional  studies  more  detailed  than  pre¬ 
viously  assumed. 

As  a  preliminary  effort  to  quantify  the  propagation  character¬ 
istics  of  seismic  waves  on  a  regional  basis,  we  have  measured  the 
attenuation  rates  of  Lg  waves  for  the  eastern  United  States,  and  the 
western  and  central  portions  of  the  Soviet  Union.  We  have  also  pro¬ 
posed  two  magnitude  scales  for  Lg  waves  and  intermediate-period  (8- 
13  sec)  Rayleigh  waves  in  the  eastern  U.S.  A  plotting  of  Mg  (from 
8-13  sec  Rayleigh  waves)  vs.  M^g  (from  0. 3-1.0  sec  Lg  waves)  for  eight 
earthquakes  and  one  underground  nuclear  explosion  in  the  eastern  U.S. 
shows  no  separation  between  the  two  populations. 

A  re-evaluation  of  the  magnitude-yield  relation  and  an  examin¬ 
ation  of  physical  parameters  which  may  be  relevant  to  the  estimated 
yield  of  underground  nuclear  explosions  were  performed.  The  prelim¬ 
inary  results  indicate  that  (i)  the  m^  vs.  yield  relation  shows  re¬ 
gional  differences  and  dependence  on  the  source  medium,  and  (ii)  the 
collapse  volume  and  the  diameter  of  the  collapsed  crater  are  usually 
proportional  to  the  estimated  yield. 


Introduction 


During  this  contract  period  we  conducted  research  on  five 
topics  which  are  directly  related  to  the  problems  of  regional  seis¬ 
mic  wave  propagation  and  earthquake-explosion  discrimination.  The 
topics  are:  (i)  a  review  of  the  available  studies  on  the  raismic 
phase  Lg,  (ii)  a  comparison  of  seismic  discrimination  methods  at 
regional  distances  in  the  U.S.  and  the  USSR,  (iii)  a  preliminary 
study  on  the  attenuation  and  magnitude-scale  of  Lg  and  intermediate- 
period  (8-13  sec)  Rayleigh  waves,  (iv)  a  preliminary  re-evaluation 
of  the  magnitude-yield  relation  and  an  examination  of  the  physical 
parameters  which  may  be  relevant  to  the  estimated  yield  of  under¬ 
ground  nuclear  explosions,  and  (v)  a  review  on  the  nature  and  re¬ 
duction  of  seismic  noise,  applicable  to  the  design  of  marine  seis¬ 
mic  systems . 

Through  the  review  on  the  properties  of  Lg  waves,  we  hope  to 
achieve  three  goals:  (i)  to  compile  and  categorize  the  available 
observations  into  accessible  format,  (ii)  to  summarize  the  theor¬ 
etical  development  in  an  overview  fashion,  and  (iii)  to  emphasize 
the  features  that  are  related  to  the  problems  of  earthquake-explosion 
discrimination.  In  this  report,  we  will  present  a  review  on  the  seis¬ 
mic  phase  Lg.  The  review  is  subdivided  into  seven  topics:  (A)  par¬ 
ticle  motion  and  dispersion,  (B)  regional  velocity,  (C)  spectral  con¬ 
tent,  (D)  wave  guide  and  mode  of  propagation,  (E)  attenuation  and 
propagation  efficiency,  (F)  magnitude-scale  based  on  Lg,  and  (G)  others 
(Sn-to-Lg  conversion,  application  to  the  earthquake-explosion  discrim¬ 
ination  problem,  and  search  for  oceanic  Lg) . 

A  comparative  study  of  regional  wave  propagation  in  the  eastern 
United  States  and  different  regions  of  the  Soviet  Union  is  presented 
in  the  second  part  of  this  Final  Technical  Report.  Four  topics  were 
selected  to  assess  the  feasibility  of  directly  comparing  the  character¬ 
istics  of  regional  seismic  waves  in  the  US  and  the  USSR,  and  to  evalu¬ 
ate  their  relative  importance  to  the  problem  of  earthquake-explosion 


discrimination.  The  topics  are:  (i)  Lg  vs.  P  amplitudes,  (ii)  Lg/P 
amplitude  ratios  as  a  function  of  distance,  (iii)  Lg  group  velocity 
at  amplitude  maxima,  and  (iv)  Lg  energy  ratios. 

To  improve  our  ability  to  discriminate  earthquakes  from  explosions 
on  a  regional  basis,  we  initiated  a  study  to  quantify  the  attenuation 
rates  of  Lg  waves  in  the  eastern  U.S.  and  the  western  and  central  por¬ 
tions  of  the  USSR.  In  addition,  we  have  suggested  a  magnitude-scale 
formula  for  Lg  waves  in  the  eastern  U.S.;  the  formula  is  very  similar 
to  the  one  proposed  by  Nuttli  (1973)  for  the  central  U.S.  Since  under¬ 
ground  nuclear  explosions  are  more  efficient  at  generating  short-period 
waves  than  long-period  waves,  we  have  begun  to  explore  the  possibility 
of  using  intermediate-period  (8-13  sec)  surface  waves  as  a  seismic  dis¬ 
criminant.  In  this  preliminary  effort,  we  have  determined  the  atten¬ 
uation  rate  and  a  magnitude  formula  appropriate  for  intermediate- 
period  Rayleigh  waves  in  the  eastern  U.S.  We  have  also  tested,  with 
negative  result,  the  potential  of  using  (from  intermediate-period 

Rayleigh  waves)  vs.  MT  as  a  discriminant.  Since  only  one  explosion 

Lg 

(SALMON)  was  used  in  the  comparison,  the  result  is  quite  inconclusive. 

In  studying  regional  seismic  wave  propagation,  we  often  encounter 
the  problem  of  how  to  calibrate  a  magnitude-yield  relation  at  regional 
distances.  This  problem,  although  quite  fundamental  in  nature,  is  by 
no  means  an  easy  one  because  a  well-determined  magnitude-yield  rela¬ 
tion  requires  a  clear  knowledge  of  (i)  the  source  size,  (ii)  the  am¬ 
plitudes  of  seismic  waves  at  different  distances,  (iii)  the  effects 
of  crustal  structure  at  the  source  and  the  receiver,  and  (iv)  the 
effects  of  the  propagation  path.  The  last  part  of  this  report  re¬ 
examines  this  relation;  it  also  describes  the  preliminary  results 
from  analyses  of  several  physical  parameters  that  are  related  to  the 
yield  of  underground  nuclear  explosions. 

A  good  understanding  on  the  nature  of  the  microseismic  noise  is 
crucial  to  the  optimization  of  signal-to-noise  ratios  for  a  marine 
seismic  system.  In  the  appendix  of  RAI's  Semi-Annual  Technical  Report, 
No.  4,  Pomeroy  (1980)  has  addressed  to  this  need  by  reviewing  the 


available  seismological  literature  on  the  subject.  In  the  review, 
he  has  focused  on  (i)  the  causal  mechanisms  and  depth-distributions 
of  propagating  and  non-propagating  seismic  noise,  (ii)  the  effects 
of  local  topography  and  geologic  structure  on  noise  amplification, 
(iii)  the  feasibility  of  installing  a  vertical  array,  and  (iv)  recom¬ 
mendations  for  the  selection  of  station  sites.  The  aforementioned 
review  is  not  reproduced  in  this  Final  Technical  Report. 
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Part  I ■  Review  of  Lq 

The  purpose  of  this  review  is  threefold:  (i)  to  provide  a 
summary  of  the  available  observations  on  Lg,  (ii)  to  present  the 
theoretical  developments  in  an  overview  fashion,  and  (iii)  to 
clarify  or  comment  on  what  appears  to  us  to  be  confusing  concern¬ 
ing  the  interpretation  of  Lg. 

The  name  Lg  was  assigned  by  Press  and  Ewing  (1952)  in  their 
pioneering  study  on  this  seismic  phase.  "L"  because  the  particle 
motion  was  predominantly  of  Love  or  transverse  type,  and  "g"  be¬ 
cause  the  wave  was  b<=  \ieved  to  propagate  in  the  granitic  layer  of 
the  crust,  and  was  t-jrefore  considered  a  surface-wave  counter¬ 
part  of  the  near-earthquake  body  waves  Pg  and  Sg.  These  authors 
summarized  the  properties  of  Lg  (for  propagation  paths  in  North 
America)  succinctly  in  the  abstract  of  their  1952  paper: 

"Surface  shear  waves  (Lg)  with  initial  period  1/2  to  6  sec¬ 
onds  with  sharp  commencements  and  amplitudes  larger  than  any  con¬ 
ventional  phase  have  been  recorded  for  continental  paths  at  dis¬ 
tances  up  to  6000  km.  These  waves  have  a  group  velocity  of  3.51 
+0.07  km/sec  and  for  distances  greater  than  20°  they  have  re¬ 
verse  dispersion.  For  distances  less  than  about  10°  the  periods 
shorten  and  Lg  merges  into  the  recognized  near-earthquake  phase  Sg." 

This  and  later  investigations  of  Lg  also  point  out  that  (i)  the 
wave  is  not  observed  after  approximately  100  km  of  propagation  in 
the  oceanic  crust,  (ii)  the  particle  motion  may  contain  a  substan¬ 
tial  amount  of  longitudinal  and  vertical  components,  and  (iii)  the 
observations  may  be  explained  by  a  collection  of  Airy  phases  of 
higher  mode  Love  and  Rayleigh  waves. 

The  terms  of  Sg  and  Lg  were  used  to  refer  to  different  waves 
in  some  earlier  studies.  Although  both  terms  referred  to  high- 
frequency  shear  waves  in  the  continental  crust,  the  distinctions 


were  based  on  differences  in  the  observed  frequency  content,  the 
distances  of  observation,  and  the  interpretation  in  their  mode  of 
propagation.  Sg,  which  is  analogous  to  its  compressional-wave 
counterpart  Pg,  referred  to  the  direct  shear  arrival  at  short  epi¬ 
center  distances;  while  Lg  referred  to  the  superposition  of  normal 
modes,  with  frequencies  slightly  lower  than  those  of  Sg,  at  epi- 
central  distances  greater  than  about  10°  (Press  and  Ewing,  1952). 
[There  has  been  considerable  confusion  concerning  the  definitions 
of  Pg  and  Sg.  These  terms  replaced  the  P  and  S  of  Mohorovicic  (1914) 
for  typographical  convenience  (page  86  of  Jeffreys,  1976)  and  the 
supposed  association  with  the  granitic  layer  of  the  crust.  While 
the  definition  of  P  referred  to  the  direct  compressional  arrival 
at  short  distances  with  a  velocity  of  about  5.5  km/sec  (cf.  Figure 
18-1  of  Richter,  1958) ,  the  original  data  was  obtained  at  distances 
over  150  km.  Explosion  data  from  California  indicated  that  direct 
compressional  arrivals  at  120  km  within  the  epicenter  had  a  veloc¬ 
ity  near  6.34  km/sec.  The  Californian  researchers  consequently 
suggested  the  notation  "p"  for  the  direct  wave  at  short  distances 
and  "P"  for  the  compressional  wave  with  a  velocity  around  5.5  km/sec 
(page  286-287  of  Richter,  1958) .  The  consensus  at  the  present  seems 
to  be  the  use  of  the  nomenclature  P  for  direct  compressional  waves 
and  the  terms  "Pn"  and  "Pg"  for  occasions  when  two  distinct  arri¬ 
vals  with  velocities  around  8. 0-8. 4  km/sec  and  5. 4-5. 7  km/sec  are 
observed.]  In  view  of  the  consensus  on  the  terminology  of  P-,  Pg-, 
and  Pn-  waves  and  the  arbitrary  distinction  between  Sg  and  Lg,  we 
are  in  favor  of  calling  the  direct  shear  arrival  "S"  and  reserving 
the  term  "Lg"  for  shear  waves  with  group  velocities  around  3.5 
km/sec  at  epicentral  distances  where  Sn  (or  the  mantle-refracted 
S)  becomes  the  first  shear  arrival.  In  this  report,  the  term  "Lg" 
will  refer  to  both  the  "Lg"  and  the  "Sg"  cited  in  earlier  seismo- 
logical  literature.  In  the  following  sections,-  we  will  attempt 
to  summarize  and  discuss  previous  studies  on  the  observations  and 
interpretations  of  the  Lg  phase.  We  have  divided  the  literature 
available  to  us  into  7  topics:  (A)  particle  motion  and  dispersion, 
(B)  regional  velocity,  (C)  spectral  content,  (D)  wave  guide  and 
mode  of  excitation,  (E)  attenuation  and  propagation  efficiency, 

(F)  magnitude-scale  based  on  Lg,  and  (G)  others. 


A.  Particle  Motion  and  Dispersion 


Press  and  Ewing  (1952)  describe  the  particle  motion  of  Lg 
in  the  following  words: 

"...During  the  first  cycles  the  waves  have  approximately 
equal  amplitudes  on  all  three  components,  but  the  transverse 
horizontal  rapidly  gains  amplitude  and  becomes  several  times 
larger  than  the  other  two  within  about  30  seconds.  Approxi¬ 
mately  one  minute  after  the  commencement  of  the  phase,  the  ampli¬ 
tude  on  the  transverse  component,  having  reached  a  value  many 
times  larger  than  that  of  S  or  SS  on  any  component,  begins  to  de¬ 
crease  gradually,  but  does  not  drop  to  a  value  comparable  with 
that  of  SS  until  about  30  minutes  later,  the  period  then  being 
of  the  order  10-14  seconds.  The  group  velocity  for  the  latter 
part  of  this  phase  is  certainly  less  than  about  2  km/sec  ,  the 
lower  limit  being  uncertain...".  As  for  Eurasian  events  recorded 
at  Uppsala  and  Kiruna,  B&th  (1954)  reports  that  the  particle  motion 
of  Lg  was  primarily  transverse  and  was  often  observed  at  two  dif¬ 
ferent  group  velocity  windows:  Lg^  at  3.54  +  0.06  km/sec  and 
Lg2  at  3.37  +  0.04  km/sec.  Lehmann  (1953)  states  that  there  was 
"considerable"  vertical  motion  involved.  All  the  authors  mentioned 
above  agreed  that  both  the  horizontal  and  the  vertical  components 
of  particle  motion  were  present  in  the  Lg  phase.  Herrin  and  Rich¬ 
mond  (1960)  used  a  ray-approach  analysis  to  explain  the  particle 
motion  of  Lg.  Their  calculations  indicate  that  a  strong  SV  type 
motion  (i.e.  with  longitudinal  and  vertical  components  of  motion) 
would  be  present  with  the  SH-type  motion  intially;  but  during  the 
later  part  of  the  wave  train  where  the  angle  of  incidence  for  the 
rays  presumably  becomes  less  steep,  energy  leakage  to  the  bottom 
layers  due  to  Sv-to-P  conversion  would  occur  and  the  SV-motion 
tends  to  decrease  faster  than  that  of  the  SH-motion.  The  results 
of  this  analysis  are  in  agreement  with  the  observations  of  Oliver 
et  al.  (1955) ,  but  do  not  agree  with  their  own  observations  at 
Dallas  for  earthquakes  in  southwestern  United  States  and  Mexico 
where  strong  SV-motion  continued  throughout  the  Lg  wave-train. 
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Herrin  and  Richmond  also  estimated  the  partitioning  of  energy  be¬ 
tween  SV  and  P  waves  at  different  angles  of  incidence;  Herrin  (1961) 
pointed  out  some  errors  in  their  partitioning  of  energy  and  cor¬ 
rected  them.  By  correlating  the  verical  component  to  the  longi¬ 
tudinal  component  of  the  Lg  particle  motion,  Sutton  et  al .  (1967) 
found  out  that  the  particle  motion  of  Lg  from  underground  nuclear 
explosions  and  small  earthquakes  tended  to  be  either  transverse 
or  mixed. 

Aside  from  the  qualitative  comparison  of  Press  and  Ewing  be¬ 
tween  the  vertical  and  horizontal  components  of  displacement,  there 
are  several  other  reports  on  their  relative  amplitudes.  For  the 
Lg  amplitudes  generated  by  the  nuclear  explosion  GNOME  in  a  salt 
mine  of  New  Mexico,  Romney  et  al.  (1962)  note  that  the  displace¬ 
ments  on  all  three  components  were  approximately  equal.  But  for 
earthquakes  in  the  northeastern  U.S. — southeastern  Canada  regions 
recorded  at  North  American  stations.  Street  (1976)  reports  that 
the  maximum  sustained  horizontal  component  of  Lg  consistently  ex¬ 
ceeded  the  vertical  component  by  a  factor  of  3.  For  all  epicentral 
distances  in  Iran,  the  resultant  horizontal  motion  of  Lg  at  1  sec 
was  usually  twice  that  of  the  vertical  component  (Nuttli,  1980a) . 

Bath  (1956) ,  however,  found  some  Lg  waves  with  no  vertical  particle 
motion  at  all. 

Although  Press  and  Ewing  (1952)  suggested  the  possibility  of 
using  higher  mode  surface  waves  to  interpret  the  Lg  phase,  Oliver 
and  Ewing  (1957)  were  the  first  to  calculate  the  dispersion  curves 
of  higher  mode  Rayleigh  waves  and  use  them  to  explain  the  longi¬ 
tudinal  and  vertical  components  of  Lg  particle  motion.  In  a  later 
paper,  Oliver  and  Ewing  (1958)  computed  the  dispersion  curves  from 
simple  earth  models  for  higher  mode  Love  and  Rayleigh  waves  and  found 
that  the  Mj-mode  (1st  shear  mode)  and  the  second  Love  mode  had  similar 
velocities  at  the  same  period,  which  may  explain  the  simultaneous  ar¬ 
rivals  of  the  vertical,  longitudinal,  and  transverse  components  of 
ground  motion  for  Lg.  Dispersion  curves  and  particle  motions  of 
higher  mode  Love  and  Rayleigh  waves  were  computed  for  realistic 
earth  models  by  Brune  and  Dorman  (1963)  ,  and  later  including  the 


effects  of  sphericity  into  the  earth  models  by  Kovach  and  Anderson 
(1964) .  Brune  and  Dorman  also  computed  synthetic  seismograms  for 
the  transverse  component  of  Lg.  The  results  of  these  authors  con¬ 
firm  the  hypothesis  of  Oliver  and  Ewing.  Knopoff  et  al .  (1973)  pre¬ 
sented  further  evidence  to  identify  the  transverse  component  of  Lg 
motion  as  higher  mode  Love  waves  by  (i)  computing  the  relative 
spectral  excitations  for  double-couple  sources  at  different  depths, 
and  { ii )  constructing  synthetic  seismograms  for  the  higher  mode 
Rayleigh  waves  and  identified  them  as  the  longitudinal  and  verti¬ 
cal  components  of  Lg  motion. 

The  particle  motion  of  the  1st  shear  mode  (M2)  was  computed 
by  Oliver  and  Ewing  (1957)  to  be  retrograde  elliptical;  the  same 
authors  later  reported  that  observations  from  an  Arctic  event 
(5/25/1950,  8:34:32;  65.5°N,  151. 5°W)  recorded  at  Palisades,  con¬ 
firmed  their  previous  theoretical  results  on  the  particle  motion 
(Oliver  and  Ewing,  1958) .  Barley  (1978)  traced  the  particle  motion 
of  higher  mode  Rayleigh  waves  (2.0  sec  _<_  T  <_  3.5  sec  )  for  the 
group  velocity  window  3.0  to  3.5  km/sec  ,  and  found  it  to  be  retro¬ 
grade  elliptical.  This  result  was  predicted  by  the  theoretical 
calculations  of  Panza  et  al.  (1972)  for  the  first  three  hiqher  Ray¬ 
leigh  modes;  these  authors  also  found  that  at  a  given  period  the 
ellipticity  (defined  as  the  ratio  of  the  longitudinal  component  of 
particle  motion  to  the  vertical  component)  increased  with  decreas¬ 
ing  mode  number.  For  a  shield  structure  with  a  low  velocity  chan¬ 
nel  (LVC)  in  the  upper  mantle,  they  found  that  at  periods  less  than 
4  sec  the  ellipticity  for  the  third  higher  Rayleigh  mode  was  greater 
or  equal  to  0.7,  whereas  the  ellipticity  for  the  fundamental  and  the 
first  two  higher  Rayleigh  modes  was  greater  or  equal  to  1.0. 

B.  Regional  Velocity 

Table  I  is  a  summary  of  Lg  velocities  which  were  published  in 
journals  and  reports  available  to  us.  Whenever  possible,  we  tried 
to  include  information  pertaining  to  the  measurements  of  the  velocity. 
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such  as  the  location  of  the  seismic  events  and  recording  stations, 
the  type  of  instrument  used  to  record  the  events  (horizontal  or 
vertical  component,  short  or  long  period,  etc.),  and  the  period 
of  the  Lg  waves  at  which  the  measurement  was  made.  Although  the 
majority  of  the  references  cited  did  not  specify  their  method  of 
measurement,  we  deduced  from  their  figures  that  most  reported 
velocities  were  measured  at  the  initial  stage  of  the  coda  when  a 
visible  change  in  wave  frequency  or  amplitude  could  be  observed, 
either  on  the  long-  or  short-period  instruments.  The  measurements 
of  Pomeroy  and  Nowak  (1978) ,  however,  were  made  at  the  amplitude 
maxima  of  the  Lg  coda  which  seemed  to  be  more  unstable.  Differences 
in  the  method  of  measurement  and  the  recording  instrument  may  ac¬ 
count  for  the  apparent  discrepancy  between  the  various  reports. 

While  measurements  at  the  beginning  of  the  coda  probably  cor¬ 
respond  to  the  Airy  phase (s)  of  higher  mode  surface  waves  with  the 
fastest  group  velocity,  measurements  at  the  amplitude  maxima  prob¬ 
ably  coincide  with  the  group  velocity  window  where  several  Airy 
phases  overlap.  Whereas  the  former  is  indicative  of  the  average 
properties  of  the  wave  guide,  the  latter  which  tends  to  be  slower 
than  the  former,  is  probably  not  only  more  diagnostic  of  the  de¬ 
tailed  structure  of  the  wave  guide  but  also  informative  concern¬ 
ing  the  relative  excitation  of  the  various  modes  at  the  source 
(Knopoff  et  al.,  1974).  We  would  like  to  explore  this  possible  as¬ 
pect  of  Lg  in  a  future  study. 

C.  Spectral  Content 

The  only  sources  known  to  us  on  the  spectral  content  of  Lg 
are  derived  from  Street  et  al.  (1976)  and  the  Soviet  seismological 
literatures  (e.g.  Antonova  et  al,  1978;  Nurmagambetov,  1974).  The 
studies  on  Lg  propagation  in  the  USSR  were  compiled  and  summarized 
in  a  report  by  Shishkevish  (1979). 

Street  et  al.  derived  their  data  from  over  300  short-period, 
vertical  component  recordings  of  78  earthquakes  in  the  central  U.S. 

*  In  the  period  range  they  analyzed  (approximately  0.05  -  10  sec  ), 
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the  amplitude  spectra  generally  indicate  a  falloff  of  (omega)  between 
the  flat  portions  at  the  long-  and  short-period  ends.  Their  spec¬ 
tra  were  corrected  for  the  effects  of  instrument  response,  but  not 
for  the  anelastic  attenuation  of  the  path. 

The  frequency  selection  seismograph  stations  (ChISS)  of  the 
USSR  have  enabled  the  spectral  analysis  of  Lg  to  become  a  routine 
procedure.  Their  results,  commonly  plotted  as  log  (A/T)  vs.  log 
( 1/T) ,  generally  display  peaks  at  short  epicentral  distances.  The 
peak  is  shifted  towards  lower  frequencies  as  epicentral  distance 
increases.  This  dependency  of  spectral  peak  on  epicentral  distance 
is  also  a  function  of  propagation  path.  In  these  studies,  the  fre¬ 
quency  ranged  from  0.3  to  approximately  20  Hz  while  the  epicentral 
distance  spanned  from  30  to  3000  km.  The  falloff  in  their  velocity 
amplitude  spectra  (i.e.  displacement  amplitude  spectra  multiplied 
by  frequency)  is  also  dependent  on  epicentral  distances:  at  epi¬ 
central  distances  around  350  km,  the  falloff  ranges  from  slightly 
greater  than  one  to  approximately  two;  whereas  at  epicentral  dis¬ 
tances  greater  than  about  1000  km,  the  falloff  remains  less  than  3. 
Since  these  measurements  of  Lg  spectral  content  did  not  take  the 
effects  of  geometrical  spreading  and  anelastic  attenuation  into 
account,  the  spectral  characteristics  measured  at  short  epicentral 
distances  were  probably  more  representative  of  the  source  spectra 
and  a  spectral  falloff  of  about  2  could  be  taken  as  representative 
of  the  source  falloff  for  the  displacement  amplitude  spectra  of  Lg 
waves.  The  high-frequency  spectral  peaks  observed  in  the  USSR  is 
probably  an  artifact  of  the  velocity  spectra  plot;  that  is,  the 
spectral  peak  will  disappear  if  the  plot  is  converted  into  a  dis¬ 
placement  amplitude  spectra. 

D.  Wave  Guide  and  Mode  of  Excitation 


Press  and  Ewing  (1952)  are,  again,  the  first  ones  to  point  out 
that  "...Lg  is  a  wave  which  is  confined  to  a  surface  or  near-surface 
layer  by  wave-guide  action..."  based  on  the  observed  velocity  and 
large  amplitudes.  Subsequent  theoretical  studies  tend  to  support 
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their  claim  although  this  conclusion  is  not  reached  without  its 
share  of  confusion.  In  a  study  of  Lg  waves  in  Eurasia,  Bath  (1954) 
observed  a  correlation  between  hypocentral  depth  and  the  energies 
contained  in  Lg^  and  I^*  That  is,  the  energy  of  Lg^  generally 
decreased  with  increasing  hypocentral  depth,  whereas  the  energy 
for  reached  a  maximum  when  the  source  depth  was  around  45  km. 

He  attributes  the  difference  in  energy  distribution  to  several 
crustal  channels  or  layers  which  transmitted  waves  at  different 
group  velocities.  This  claim,  although  sound  when  interpreted 
in  terms  of  Airy  phases  with  different  group  velocities,  led  to 
two  unexpected  results  when  viewed  from  the  perspective  of  chan¬ 
nel  waves.  Firstly,  terminologies  for  waves  which  supposedly 
propagated  in  different  channels  of  the  crust  and  upper  mantle 
proliferated  (e.g.  B&th,  1958).  Secondly,  several  low-velocity 
channels  in  the  crust  and  upper  mantle  came  to  be  used  as  ex¬ 
planations  for  the  efficient  propagation  of  the  various  channel 
waves  (Gutenberg,  1955;  B&th,  1956,  1958) . 

Based  on  the  dispersion  curves  of  higher  mode  Love  and  Ray¬ 
leigh  waves,  Oliver  and  Ewing  (1957,  1958),  Brune  and  Dorman  (1963), 
and  Kovach  and  Anderson  (1964)  found  it  possible  to  explain  the 
frequency  content  and  the  group  velocity  of  Lg  waves  by  using  the 
Airy  phases  of  the  higher  modes.  Kovach  and  Anderson  (1964)  also 
point  out  that  the  modes  observed  "...depend  on  the  period  range 
being  studied  and  the  depth  of  the  source..."  and  that  variations 
in  the  velocity  and  period  of  the  observed  Lg  depended  on  the  po¬ 
sitions  of  the  Airy  phase,  which  in  turn  depended  on  the  elastic 
parameters  of  the  propagation  path.  If  the  interpretation  of  Lg 
waves  as  superpositions  of  higher  mode  surface  waves  is  correct, 
then  we  would  expect  an  additional  dependence  on  the  source  radi¬ 
ation  pattern.  At  periods  greater  or  equal  to  5  sec  ,  radiation 
patterns  of  the  first  higher  Love  and  Rayleigh  modes  compare 
favorable  with  calculated  results  (Mitchell,  1973,  a,b) .  The  ob¬ 
servations  of  Sutton  et  al.  (1967)  on  short-period  (0. 5-2.0  sec  ) 

Lg  waves,  however,  indicate  that  "...there  seems  to  be  no  system¬ 
atic  difference  in  the  short-period  energy  radiation  pattern  be- 
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tween  the  underground  nuclear  explosions  and  the  earthquakes..." 
and  that  the  pattern  of  the  energy-contours  (or  contours  based 
on  the  maximum  amplitude)  could  be  better  explained  by  a  cor¬ 
relation  with  the  major  tectonic  provinces  of  the  United  States. 
Since  the  modal  composition  of  Lg  at  short  periods  is  a  combin¬ 
ation  of  many  higher  modes,  the  observed  amplitudes  may  not  be 
diagnostic  of  the  radiation  pattern  of  the  individual  modes. 

Also,  scattering  is  probably  more  important  for  short-period 
waves  and  its  effects  more  likely  to  mask  any  azimuthal  pattern 
that  may  be  present. 

Panza  et  al.  (1972)  showed  that  the  collection  of  higher 
mode  Rayleigh  waves  could  be  separated  into  a  family  of  crustal 
waves  and  a  family  of  channel  waves  in  a  structure  containing 
even  a  slight  low-velocity  channel  (LVC)  in  the  upper  mantle. 

As  it  is  implied  by  the  name,  channel  waves  have  most  of  the 
energy  in  the  LVC  and  have  essentially  zero  energy  at  the  sur¬ 
face.  Crustal  waves,  on  the  other  hand,  have  most  of  their 
energy  in  the  crust;  consequently,  only  the  fundamental  mode 
and  the  crustal  waves  need  to  be  considered  for  the  excitation 
of  Rayleigh  waves.  Knopoff  et  al.  (1973)  demonstrated  that  higher 
mode  Love  waves  could  similarly  be  divided  into  crustal  waves 
and  channel  waves.  For  a  structure  without  any  LVC,  the  whole 
suite  of  higher  mode  Love  and  Rayleigh  waves  has  to  be  taken 
into  account  for  the  ground  motion  of  the  Lg  waves. 

Knopoff  et  al.  (1974)  further  establish  that  the  group  ve¬ 
locity  and  the  periods  of  the  Lg  stationary  phase  could  be  diag¬ 
nostic  for  the  crustal  thickness  and  the  shear  velocity  in  the 
crust  and  the  upper  mantle.  In  general,  as  the  crustal  thickness 
increased,  both  the  group  velocity  of  the  late-arriving  Lg 

stationary  phases,  U  .  ,  and  the  period  at  U  .  ,  T  .  ,  tended  to 

min  mm  mm 

increase.  Increasing  the  crustal  velocity  while  keeping  all  other 

parameters  constant  would  tend  to  decrease  T  .  ,  but  increase 

mm 

^min'  t*ie  magnitude  of  Lg-excitation,  and  the  general  period-content 
of  the  Lg  waves.  These  authors  also  demonstrate  that  (i)  for  thick- 


nesses  of  the  upper  mantle  lid  greater  than  20-25  km,  Lg  is  in¬ 
sensitive  to  changes  in  its  thickness,  and  (ii)  Lg  is  insensi¬ 
tive  to  the  velocity  in  the  upper  mantle  LVC.  Panza  and  Calcan- 
nile  (1975)  point  out  that  higher  mode  contribution  becomes  more 
significant  as  the  period  decreases  and/or  as  the  hypocentral 
depth  increases. 

As  for  the  low-velocity  channel  in  the  crust  and/or  upper 
mantle,  Oliver  and  Ewing  (1958)  concluded  that  it  was  not  neces¬ 
sary  to  explain  the  characteristics  of  the  Lg  phase.  Knopoff 
et  al.  (1973)  and  Panza  and  Calcaqnile  (1975) ,  based  on  more  nodes  ex 
tending  to  shorter  periods,  reached  the  same  conclusion  concern¬ 
ing  the  Love-  and  Rayleigh-type  motions  of  the  Lg  phase,  re¬ 
spectively. 

Most  of  the  investigators  mentioned  in  this  section  would 
probably  maintain  that  the  characteristics  of  Lg  can  be  explained 
by  the  anelastic  attenuation  of  the  crust-mantle  layers,  the  fre¬ 
quency  response  of  the  seismograph  system,  and  the  superposition 
of  higher  mode  surface  waves.  Ruzaikin  et  al.  (1977)  ,  on  the  other 
hand,  state  that  they  "...remain  unconvinced  that  normal  modes 
will  allow  useful  interpretation  of  Lg  when  more  detailed  data  on 
its  structure  are  obtained..."  and  suggest  that  lateral  hetero¬ 
geneity  had  a  key  role  in  shaping  the  characteristics  of  the  ob¬ 
served  Lg.  Their  argument  was  based  on  the  discrepancy  between 
calculations  from  higher  mode  surface  waves  which  predicted  the 
duration  of  Lg  to  be  confined  in  the  group  velocity  windows  of 
approximately  3. 5-3.1  km/sec.,  and  observations  of  the  Lg  phase 
which  indicated  that  its  amplitude  was  significant  in  the  group 
velocity  window  3. 5-2. 8  km/sec.  Oceanic  Rayleigh  waves  of  the 
fundamental  mode  (T  12  sec  )  also  exhibit  similar  "stretching" 
in  duration.  These  waves  have  nevertheless  been  instrumental  in 
shaping  our  present  understanding  concerning  the  oceanic  structure. 
Thus,  while  we  share  the  belief  with  Ruzaikin  et  al.  that  hetero¬ 
geneities  in  the  propagation  path  are  important  in  shaping  the 
waveform  of  Lg,  we  also  believe  that  the  normal  mode  theory,  when 
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supplemented  with  theories  or  methods  which  can  take  heterogeneity 
in  the  path  into  consideration  (e.g.  the  scattering  theory  of  Aki, 
1969) ,  will  serve  to  improve  the  explanation  for  the  Lg  phase. 


E.  Attenuation  and  Propagation  Efficiency 

This  section  deals  with  the  measurement  of  amplitude-diminution 
as  a  function  of  epicentral  distance;  the  title  of  the  section  re¬ 
flects,  respectively,  the  quantitative  and  qualitative  aspects  of 
it.  The  former  refers  to  the  rate  of  anelastic  absorption  of  the 
wave's  kinetic  energy  per  unit  distance,  while  the  latter  provides 
a  descriptive  measure  for  the  efficiency  of  the  medium  in  trans¬ 
mitting  Lg  waves. 

In  seismological  literature,  attenuation  is  usually  measured 
in  terms  of  the  attenuation  coefficient,  gr  ,  or  the  attenuation 
quality  factor,  Q.  These  two  quantities  can  be  related  via  the 
following  equation: 

*  -  “q"d~  <» 

where  f  and  u  are  the  frequency  and  the  velocity  of  the  wave,  re¬ 
spectively.  For  Lg  waves,  measurements  of  'Q'  and  Q,  compiled  in 
Table  II,  have  been  obtained  by  three  approaches:  (i)  time-domain, 
(ii)  frequency-domain,  and  (iii)  coda. 

The  time  domain  approach  entails  three  steps:  (i)  measure 
the  wave  amplitude  at  different  epicentral  distances,  (ii)  correct 
the  amplitudes  for  the  effect  of  geometrical  spreading,  and  (iii) 
estimate  the  ^  or  Q  that  would  explain  the  falloff  of  the  ampli¬ 
tude  in  relation  to  distance.  Nuttli  (1975,  1978,  1980  a,b)  and 
Street  (1976)  chose  to  combine  steps  (iii)  and  (ii)  together,  and 
compared  the  observed  amplitudes  directly  with  curves  that  include 
the  effects  of  geometrical  spreading  and  different  degrees  of  atten¬ 
uation.  The  frequency-domain  approach  has  the  advantage  of  being 
able  to  take  the  source  radiation  pattern  into  account.  The  pro¬ 
cedure  used  by  Mitchell  and  coworkers,  who  have  been  the  primary 
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advocates  of  this  approach  on  higher  mode  surface  waves,  is  simi¬ 
lar  to  that  employed  for  the  study  of  the  fundamental  mode  (Tsai 
and  Aki,  1969).  Again,  three  steps  are  involved  in  this  procedure: 
(i)  determine  the  amplitude  spectra  for  the  fundamental  and  higher 
mode  surface  waves  by  applying  a  f requency-velocity  filter  (e.g. 
the  multiple-filter  technique  of  Dziewonski  et  al.,  1969)  (ii) 
calculate  the  attenuation  coefficient  that  would  produce  the  best 
fit  between  the  observed  amplitudes  and  the  radiation  pattern  com¬ 
puted  at  each  period.  To  date,  this  approach  has  been  limited  to 
the  analysis  of  the  fundamental  and  the  1st  higher  mode  (Mitchell, 
1973  a,b;  Cheng  and  Mitchell,  1980) .  The  coda  approach,  which  was 
derived  from  the  scattering  theory  of  surface  waves  (Aki,  1969) , 
has  been  applied  successfully  to  data  from  narrow-band  seismographs 
to  establish  (i)  scaling  laws  for  local  earthquakes,  and  (ii)  esti¬ 
mates  of  regional  Q  (Aki  and  Chouet,  1975;  Chouet  et  al.,  1978; 
Rautian  and  Khalturin,  1978)  .  Herrmann  and  coworkers  recently 
modified  this  method  for  data  derived  from  broadband  seismographs. 
They  estimated  the  regional  Q  from  Lg  waves  by  measuring  (i)  the 
predominant  frequency  in  the  coda  as  a  function  of  time,  and  (ii) 
the  coda  shape  (Herrmann,  1980;  Singh  and  Herrmann,  1979) . 

The  propagation  efficiency  of  a  region  is  usually  estimated 
by  measuring  the  frequency  content  and  wave  amplitude  (usually  in 
relation  to  the  level  of  the  ambient  noise  or  the  amplitude  of  an¬ 
other  phase);  in  general,  three  terms:  clear,  weak,  and  none  are 
used  to  describe  the  amplitude  of  the  Lg  phase.  "Clear"  usually 
refers  to  an  impulsive,  large-amplitude,  high-frequency  arrival; 
"weak"  refers  to  a  drawn-out,  small,  low-frequency  arrival;  and 
"none"  is  indicative  of  completely  inefficient  Lg  propagation. 
Although  different  authors  have  set  their  standards  for  clear  and 
weak  Lg  somewhat  differently,  their  conclusions  concerning  the 
propagation  efficiency  of  a  given  region  are,  surprisingly,  quite 
uniform.  A  list  of  regional  studies  on  the  propagation  efficiency 
of  Lg  is  compiled  in  Table  III. 
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In  interpreting  the  inefficient  propagation  of  Lg  in  the 
Tibetan  plateau,  Ru2aikin  et  al .  (1977)  proposed  two  explanations 
which  are  probably  applicable  to  most  areas  with  major  tectonic 
boundaries.  Firstly,  a  disruption,  termination,  or  vertical  dis¬ 
placement  of  wave  guide  (which  is  either  the  entire  crust  or  part 
of  it)  will  seriously  affect  the  propagation  efficiency  of  Lg  waves; 
secondly,  high  attenuation  in  the  crust  will  also  be  able  to  affect 
the  ability  to  transmit  Lg.  The  ocean-continent  boundary  is  prob¬ 
ably  a  disruption  or  termination  of  the  wave  guide  for  Lg;  disap¬ 
pearance  of  Lg  waves  after  crossing  approximately  100  km  of  oceanic 
structure  is  a  well  documented  observation  (e.g.  Press  and  Ewing, 
1952;  Bath,  1954;  etc.).  This  peculiar  property  of  Lg  waves  to 
propagate  only  in  the  continental  crust  was  used  by  Oliver  et  al. 
(1955)  to  map  the  continental  structure  in  the  Arctic  regions. 

B4th  (1956)  and  Gutenberg  (1955)  report  that  the  Lg  phase 
was  weakened  or  disappeared  when  crossing  recent  mountain  chains. 
Shishkevish  (1979) ,  in  his  compilation  of  studies  on  Lg  propaga¬ 
tion  in  the  Soviet  Union,  also  notes  that  the  Lg  phase  was  atten¬ 
uated  when  crossing  Tien  Shan,  Pamir-Hindu  Kush,  and  the  Himalayas. 
He  also  point  out  that  "...the  propagation  of  Lg  across  the  Tien 
Shan  is  less  efficient  when  paths  are  more  oblique  to  the  trend 
of  the  range  than  when  they  are  perpendicular  to  it...".  Uni¬ 
formity  of  the  structure  (Chinn  et  al.,  1980)  and  the  complexity 
of  geology  (Street,  1976)  in  the  propagation  path  are  also  con¬ 
sidered  important  in  determining  the  attenuation  of  the  Lg  am¬ 
plitude.  In  summary,  the  presence  of  a  uniform,  high-Q  wave 
guide  is  essential  for  the  efficient  propagation  of  Lg;  in  the 
case  of  a  non-uniform  or  low-Q  wave  guide,  the  degree  of  non¬ 
uniformity  of  the  wave  guide  and  the  length  of  propagation  in  it 
are  both  important  in  determining  the  fraction  of  Lg-energy  that 
will  be  observed. 

F.  Magnitude-Scale  Based  on  Lg 


Since  Lg  is  often  found  to  be  the  largest  phase  at  regional 


distances,  it  is  natural  that  a  magnitude-scale  based  on  Lg  am¬ 
plitude  would  become  important  to  studies  on  regional  seismicity. 
Based  on  LRSM  reports  from  78  underground  nuclear  explosions. 

Baker  (1970)  proposed  a  general  formula  of  the  form, 

MLg  =  log10  (A/T)  +  Q(T,a)  +  S  (T)  (2) 

to  calculate1  the  magnitude-scale  from  Lg  amplitudes.  Q(T,A) 
represents  a  correction  term  for  the  attenuation,  and  S (T)  is  a 
term  for  station  correction.  Baker  obtained  an  expression  for 
Q(T,A),  as  a  sixth  degree  polynomial  of  distance,  by  minimizing 
the  difference  between  log^g(A/T)  and  the  reported  m^  for  each 
event;  he  also  assigned  tentative  corrections  for  each  station, 
calculated  by  Baker  indicates  less  scatter  than  the  reported 

V 

Nuttli  (1973)  formulated  a  magnitude  scale  for  Lg  while 
studying  its  attenuation  in  the  eastern  United  States.  He  assumed 
that  the  term  Q  (T,  A )  in  equation  (2)  has  the  form  C(T,A)  log^QA, 
and  subsequently  found  two  magnitude  formulae,  applicable  at  dif¬ 
ferent  distance  ranges,  for  1-sec  Lg  of  "sustained"  (3  or  more 
cycles)  amplitudes. 

MLg  =  3,75  +  0,9  log10A+  logl0  (A/T)  0.5o1a1  4° 

=  3.30  +  1.66  log1()A  +  log1Q  (A/T)  4°^-^^-30° 

Street  (1976)  and  Bollinger  (1979),  respectively,  found  Nuttli's 
formulae  to  be  applicable  in  northeastern  and  southeastern  North 
America,  provided  that  the  maximum  distance  is  limited  to  approx¬ 
imately  2000  km. 


Street  et  al.  (1976),  on  the  other  hand,  assumed  C(T,A)  to  be 
known  and  then  specified  S (T)  such  that  the  magnitude  scales  at  dif¬ 
ferent  periods  were  set  equal  for  an  m^  =  1.5  event.  For  an  m^  =  2.5 
event,  the  magnitude  calculated  at  0.1  sec.  according  to  their  formu¬ 
lation  would  be  1.8,  and  the  discrepancy  between  m^  and  m^  ^  increased 
rapidly  with  increasing  m^.  Since  there  is  no  implicit  or  explicit 
reasoning  behind  the  assumption  of  a  known  C(T, A) ,  we  are  inclined 
towards  the  procedure  of  determining  C(T,A)  experimentally  and  then 
calculating  the  S(T)  so  that  a  uniform  magnitude  would  be  obtained 
at  all  periods. 
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G.  Others 


Sn  to  Lg  conversion  appears  to  occur  near  the  margin  of  the 
American  continents.  For  events  from  the  West  Indies  and  Mexico 
recorded  at  North  American  stations,  Isacks  and  Stephens  (1975) 
identified  the  prominent  phases  which  arrived  after  Sn  as  pos¬ 
sibly  a  converted  Lg  at  the  continental  margin.  Chinn  et  al. 

(1980)  observed  similar  conversions  for  events  in  the  Nazca 
Plate  recorded  at  South  American  stations.  In  neither  of  the 
studies  was  any  Lg  to  Sn  conversion  observed. 

A  number  of  investigators  have  explored  the  possibility  of 
using  the  ratio  of  Lg-amplitude  to  P-amplitude  as  a  discriminant 
for  the  earthquake  and  the  underground  explosion  populations. 

This  possibility  was  tested  by  Pomeroy  and  Nowak  (1979) ,  Pomeroy 
(1980),  Nuttli  (1980  b)  ,  and  Gupta  et  al.  (1980)  for  propagation 
paths  in  western  and  central  Soviet  Union,  and  by  Pomeroy  and 
Nowak  (1979)  and  Pomeroy  (1980)  for  propagation  paths  in  eastern 
and  western  United  States,  respectively.  Their  findings  indi¬ 
cate  a  tendency  for  the  Lg  to  P  amplitude  ratios  to  be  greater 
than  1.0  for  earthquakes  and  less  than  1.0  for  underground 
nuclear  explosions.  The  ratios,  however,  appear  to  be  strongly 
dependent  on  the  epicentral  distance  and  the  regional  attenua¬ 
tion  in  the  propagation  paths  and  therefore  cannot  be  used  re¬ 
liably  as  a  discriminant  between  explosions  and  earthquakes. 

Contrary  to  higher-mode  surface  waves  in  continental  struc¬ 
tures,  higher-mode  Love  waves  in  sediment-covered  oceanic  struc¬ 
tures  do  not  form  a  coherent  family  of  arrivals  at  short  periods 
(Knopoff  et  al.,  1979).  This  phenomenon  can  serve  to  explain  the 
absence  of  Lg  waves  in  the  oceanic  structure.  These  authors  also 
point  out  that  since  a  large  fraction  of  the  shear  energy  at  the 
stationary  phases  of  higher-mode  Love  waves  is  concentrated  in 
the  sedimentary  layer,  absorption  by  the  low-rigidity  sediment  and 
scattering  due  to  variations  in  its  thickness  can  account  for  the 
rapid  attenuation  of  the  higher-mode  Love  waves  in  oceanic  structures. 
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Part  II.  Seismic  Discrimination  Methods  at  Regional  Distances 

A.  Propagation  Characteristics 

For  propagation  paths  within  eastern  North  America  (ENA)  Lg 
is  commonly  the  phase  with  the  largest  amplitudes  on  conventional 
short-period  seismograms  (with  peak  response  at  less  or  near  1  Hz) . 
Propagating  at  a  group  velocity  of  approximately  3.5  km/sec.,  the 
recorded  Lg  has  (i)  predominant  frequencies  of  1  to  3  Hz,  (ii) 
particle  motion  in  all  three  components:  transverse,  longitudi¬ 
nal,  and  vertical,  and  (iii)  maximum  amplitudes  six  to  ten  times 
larger  than  those  of  P  waves.  Works  of  Ruzaikin  et  al.  (1977) 
and  Antonova  et  al.  (1978)  indicate  that  Lg  propagation  in  the 
eastern  USSR  seems  comparably  efficient,  with  Lg  amplitudes  sig¬ 
nificantly  larger  than  the  P  amplitudes.  Furthermore,  Pg  propa¬ 
gation  is  generally  inefficient  in  both  regions.  In  contrast, 
for  propagation  paths  within  western  North  America,  as  well  as  in 
the  western  and  central  portions  of  the  USSR  Lg  is  observed  to 
have  roughly  equal  amplitudes  as  P,  but  Pg  appears  relatively 
more  prominent  at  regional  distances.  These  observations: 
large  Lg  and  small  Pg  or  vice  versa,  seem  to  imply  a  relation 
between  Lg-  and,  Pg-  propagation  and  the  crustal  structures  along 
the  propagation  paths.  Following  Haskell's  (1966)  interpretation 
that  (i)  the  attenuation  of  short-period,  continental  crustal  P 
waves  may,  to  a  large  degree,  be  explained  as  leakage  of  energy 
to  the  layers  beneath  the  waveguide,  and  that  (ii)  low  leakage 
can  generally  be  associated  with  low  near-surface  velocities,  we 
are  in  favor  of  explaining  the  relative  amplitudes  of  Pg  in  terms 
of  velocity  contrasts  at  the  lower  interface  of  the  Pg  waveguide 
(presumably  the  Moho)  and/or  near-surface  structures  in  the  propa¬ 
gation  paths.  The  amplitudes  of  Lg  may  also  be  related  to  these 
velocity  structures  in  the  waveguide,  but  the  exact  relation  is 
not  clear. 

B.  Lg-  Amplitude  vs.  P-  Amplitude 

In  this  section,  we  present  the  quantitative  relation  between 
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the  maximum  amplitudes  of  P-  and  Lg-  waves.  Measurements  of  wave 
amplitudes  were  made  from  seismograms  reocrded  by  the  short-period 
instruments  of  World-Wide  Standard  Seismograph  Network  (WWSSN) , 
Northeastern  U.S.  Seismic  Network  (NEUSSN)  operated  by  Lamont- 
Doherty  Geological  Observatory,  and  Long  Range  Seismic  Measure¬ 
ments  (LRSM) .  The  response  curves  for  these  instruments  are  shown 
in  Figure  1.  For  the  earthquakes  and  explosions  listed  in  Table  IV 
and  V  and  plotted  in  Figures  2  to  3,  measurements  for  the  maxi¬ 
mum  P-  and  Lg-  wave  amplitudes  were  made  at  the  stations  where 
both  waves  could  be  identified.  The  ground  motions  were  then 
calculated  from  the  measured  amplitudes  by  correcting  for  the  in¬ 
strument  magnification.  Comparisons  of  ground  motions  for  P-  and 
Lg-  waves  are  presented  in  Figures  4  to  8  for  the  eastern  U.S., 
eastern  USSR,  western  U.S.,  central  and  western  USSR,  respectively. 
The  results  for  the  eastern  USSR  were  taken  from  Ruzaikin  et  al. 
(1977)  and  Antonova  et  al.  (1978).  Since  no  scale  was  given  for 
the  seismograms  reproduced  in  these  two  papers,  the  wave  amplitudes 
in  Figure  5  represents  the  record  amplitude  in  millimeters. 

Before  we  proceed  to  discuss  the  results,  we  would  like  to 
point  out  that  for  events  located  in  western  and  central  USSR  the 
source  and  the  receiver  are,  in  general,  separated  by  tectonic 
boundaries,  whereas  for  eastern  and  western  U.S.  and  eastern  USSR 
the  observations  were  usually  made  within  the  same  tectonic  prov¬ 
ince  as  the  sources.  The  inclusion  of  tectonic  boundaries  in  the 
propagation  path  may  introduce  significant  effects  on  both  the 
amplitude  and  the  phase,  or  in  short,  the  waveform  of  the  seismic 
phase;  consequently,  we  believe  that  direct  comparisons  between 
results  from  the  different  regions  should  be  made  with  great  caution. 

An  examination  of  Figures  4  through  8  indicates  that  (i)  Lg- 
amplitudes  are  much  larger  than  P-  amplitudes  in  the  eastern  parts 
of  the  United  States  and  the  Soviet  Union,  (ii)  for  propagation 
paths  in  western  U.S.,  western  and  central  USSR  Lg-  and  P-  ampli¬ 
tudes  are  comparable,  and  (iii)  the  amplitude  ratios  of  Lg  to  P 
for  earthquakes  are  somewhat  larger  than  those  for  explosions. 
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Figure  3.  Location  of  USSR  events 
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Figure  4.  Lg-  amplitude  vs.  P-  amplitude  in  the  eastern  U.S 
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Lg-amplitude  vs.  P- amplitude  in  the  eastern  USSR. 
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Figure  6.'  La-  amplitude  vs.  P-  amplitude  in  the  western  U.S 
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Lg-  amplitude  vs,  p-  amplitude  in  the  central  USSR 
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Figure  8.  Lg-  amplitude  vs.  P-  amplitude  in  the  western  USSR. 
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Table  IV 


Earthquakes  Used  in  This  Study 
Eastern  U.S. 


Date 

Origin  Time 

Location 

Latitude 

Longitude 

Comments 

06/15/73 

01:09:05 

Maine-NH 

Quebec 

Border 

45.390° 

71.000° 

5.2 

%=4.9 

01/08/74 

01:12:37.4 

Tennessee 

36.20° 

89.39° 

4.1 

4 . 3  (S ) 

02/15/74 

22:35:44.7 

Arkansas 

34.05° 

93.13° 

4.2 

3.6 (S) 

04/03/74 

23:05:02.5 

S. Illinois 

38.59° 

88.09° 

4.5 

4.7  (S) 

06/05/74 

08:06:11.3 

S. Illinois 

38.62° 

89.94° 

4.0 

3 . 6  (S) 

06/13/75 

22:40:27.2 

Missouri 

36.54° 

89.68° 

4.3 

07/09/75 

14:54:15.1 

Minnesota 

45.67° 

96.04° 

4.6 

4 . 3 (S) 

07/12/75 

12:37:16 

Maniwaki 

46.467° 

76.222° 

4  ’  lmN 

08/29/75 

04:22:51.9 

Alabama 

33.82° 

86.60° 

3.5 

4 . 4 (S) 

10/23/75 

21:17:48.2 

Manicouagan 

49.689° 

68.822° 

4 "  Oh'n 

10/23/76 

20:58:18 

St.  Simeon/ 
Quebec 

47.492° 

69.474° 

4  *  4mN 

10/15/63 

12:28:58.4 

Southern 

Quebec 

46.6° 

77.6° 

<3 

10/16/63 

15:31:01.8 

Southern 

New  England 

42.5° 

70.8° 

<3 

10/10/63 

14:59:52.5 

Virginia 

39.8° 

78.2° 

<3 

05/04/63 

21:01:35.9 

S . Carolina 

32.3° 

79.7° 

<3 

12/04/63 

21:32:35.1 

Northern 

New  England 

43.6° 

71.5° 

<3 

12/05/63 

06:51:02.5 

Kentucky 

37.2° 

87.0° 

<3 

02/18/78 

14:48:25.3 

Canada 

46.31° 

74.37° 

4.2 

08/14/65 

13:13:56.6 

S. Illinois 

37.23° 

89.28° 

3.8 

Western  USSR 

03/02/66 

02:37:03 

Caucasus 

Mtns. 

43.03° 

45.71° 

4.9 

02/21/70 

07:09:15 

Urals 

59.40° 

59.80° 

4.0+.5 

03/21/76 

22:39:40.2 

Central 

Kazakh 

42.97° 

69.89° 

3.9 

04/02/76 

17:52:28.3 

E. Caucasus 

42.99° 

45.09° 

4.5 

Table  iv 


Earthquakes  Used  in  This  Study 
Western  USSR 


Date 

Origin  Time 

Location 

Latitude 

Longitude 

% 

Comments 

04/08/76 

22:54:17. 8 

Uzbek 

40.487° 

63.650° 

4.7 

04/29/76 

23:23:15.7 

Turkey-USSR 

Border 

40.977° 

42.874° 

4.8 

10/25/76 

08:39:46.4 

Europe-USSR 

Border 

59.157° 

23.725° 

4.5 

Central  USSR 


04/08/76 

22:54:17.8 

Gazli 

40.487° 

63.650° 

04/12/76 

16:12:58.9 

Gazli 

40.456° 

63.610° 

04/17/76 

20:21:47.2 

Gazli 

40.446° 

63.686° 

04/18/76 

22:37:39.7 

Gazli 

40.265° 

63.812° 

04/21/76 

22:33:29.8 

Gazli 

40.550° 

63.846° 

Eastern  USSR 

05/22/73 

02:15:04 

52.9° 

89.5° 

02/27/72 

22:15:03 

55.1° 

93.1° 

04/30/71 

15:45:12 

46.4° 

96.6° 

03/25/72 

05:58:05 

44.9° 

101.0° 

02/04/72 

03:34:48 

53.1° 

107.8° 

12/18/71 

22:23:48 

56.6° 

114.0° 

01/15/72 

18:08:04 

58.2° 

120.7° 

11/25/72 

13:42:34 

56.3° 

123.6° 

08/09/72 

20:51:50 

56.9° 

127.7° 

12/29/73 

14:41:31 

44.7° 

82.8° 

07/07/73 

11:41:25 

40.0° 

77.4° 

03/15/73 

23:24:25 

37.4° 

77.8° 

02/23/73 

10:45:08 

37.9° 

86.9° 

07/16/73 

19:45:43 

35.  3° 

86.4° 

10/13/74 

21:29:47 

34.8° 

87.4° 

12/30/72 

23:54:09 

34.0° 

87.6° 

02/04/72 

14:08:20 

30.6° 

84.4° 

08/10/72 

21:06:41 

32.5° 

93.7° 

07/17/71 

15:00:53 

26.2° 

93.3° 

I 

t 
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Table  V 


Explosions  Used  in  This  Stuch 
Western  USSR 


Estimated 


Date 

Oriqin  Time 

Location 

Latitude 

Longitude 

Yield  i 

03/23/71 

06:59:56 

Urals 

61.29° 

56.47° 

51 

07/02/71 

17:00:02 

Urals 

67.66° 

62.00° 

7 

07/10/71 

16:59:59 

Urals 

64.17° 

55.18° 

27 

09/19/71 

11:00:07 

Urals 

57.78° 

41.10° 

4 

10/04/71 

10:00:02 

Urals 

61.61° 

47.12° 

11 

10/22/71 

05:00:00 

Urals 

51.57° 

54.54° 

34 

12/30/71 

06:20:58 

Semi- 

palatinsk 

49.75° 

78.13° 

“  — 

07/09/72 

06:59:58 

N. Black  Sea 

49.78° 

35.04° 

6 

08/20/72 

02:59:58 

N. Caspian 
Sea 

49.46° 

48.18° 

87 

09/21/72 

09:00:01 

N. Caspian 
Sea 

52.13° 

51.99° 

21 

10/03/72 

08:59:58 

N.W. Caspian 
Sea 

46.85° 

45.01° 

88 

11/24/72 

09:59:58 

W. Kazakh¬ 
stan 

51.84° 

64.15° 

20 

09/30/73 

04:59:57 

Urals 

51.61° 

54.58° 

22 

10/26/73 

05:59:58 

Urals 

53.66° 

55.38° 

7 

07/08/74 

06:00:02 

Urals 

53.80° 

55.20° 

— 

08/29/74 

15:00:00 

Urals 

67.23° 

62.12° 

20 

06/09/76 

03:02:58 

E. Kazakh 

50.02° 

79.08° 

— 

07/04/76 

02:56:58 

E. Kazakh 

49.91° 

78.95° 

— 

07/29/76 

02:57:00 

E. Kazakh 

50.00° 

78.00° 

— 

03/29/77 

03:56:58 

E. Kazakh 

49.79° 

78.15° 

— 

Central  USSR 

09/29/68 

03:42:58 

E. Kazakh 

49.77° 

78.19° 

11/09/68 

02:53:58 

E. Kazakh 

49.79° 

78.04° 

12/18/68 

05:01:57 

E. Kazakh 

49.72° 

78.06° 

03/07/69 

08:26:58 

E. Kazakh 

49.81° 

78.15° 

07/04/69 

02:46:57 

E. Kazakh 

49.75° 

78.19° 

*  9 
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Table  v 


Explosions  Used  in  This  Study 


Central  USSR 


Date 

Oriqin  Time 

09/11/69 

04:01:57 

10/01/69 

04:02:58 

07/21/70 

03:02:57 

11/04/70 

06:02:57 

12/17/70 

07:00:57 

04/25/71 

03:32:58 

05/25/71 

04:02:58 

10/09/71 

06:02:57 

02/10/72 

05:02:57 

03/28/7 

04:21:57 

11/02/72 

01:26:58 

02/16/73 

05:02:58 

07/10/73 

01:26:58 

07/23/73 

01:22:58 

05/31/74 

03:26:57 

12/27/74 

05:46:57 

02/02/75 

05:32:58 

06/08/75 

03:26:58 

12/13/75 

04:56:57 

12/25/75 

05:16:57 

12/06/69 

07:02:57 

12/12/70 

07:00:57 

12/23/70 

07:00:57 

04/11/72 

06:00:05 

Location 

Latitude 

E. Kazakh 

49.70° 

E. Kazakh 

49.81° 

E. Kazakh 

49.95° 

E. Kazakh 

49.97° 

E. Kazakh 

49.73° 

E. Kazakh 

49.82° 

E. Kazakh 

49.80° 

E. Kazakh 

50.00° 

E. Kazakh 

49.99° 

E. Kazakh 

49.73° 

E. Kazakh 

49.91° 

E. Kazakh 

49.83° 

E. Kazakh 

49.78° 

E. Kazakh 

49.99° 

E. Kazakh 

49.95° 

E. Kazakh 

49.96° 

E. Kazakh 

49.82° 

E. Kazakh 

49.76° 

E. Kazakh 

49.80° 

E. Kazakh 

50.04° 

E. Kazakh 

43.83° 

E. Kazakh 

43.85° 

E. Kazakh 

43.84° 

E. Kazakh 

37.37° 

Estimated 
Longitude  Yield  (kt) 

78.11° 

78.21° 

77.75° 

77.79° 

78.13° 

78.07° 

78.21° 

77.70° 

78.89“ 

78.19° 

78.84° 

78.23° 

78.06° 

78.85° 

78.84° 

79.05° 

78.08° 

78.09° 

78.20° 

78.90° 

54.78“ 

54.77° 

54.85° 

62.00° 


Comments 


♦Information  on  the  earthquakes  and  explosions  in  the  western  North  America 
is  available  on  request  from  RAI * 


Significant  overlap  between  the  two  populations,  however,  pre¬ 
vents  the  amplitude-ratio  method  from  becoming  an  effective  dis¬ 
criminant.  The  last  observation  is  in  agreement  with  Nuttli 
(1980  b)  for  events  in  western  and  central  Asia,  but  contrary 
to  the  conclusion  of  Gupta  et  al.  (1980)  for  propagation  paths 
in  western  Russia. 

Additional  results  on  the  amplitude  ratios  of  Lg  to  P  waves 
as  a  function  of  distance  are  presented  in  Figures  9a  to  9e. 

These  plots  were  taken  from  the  published  results  of  Soviet  in¬ 
vestigators  as  compiled  by  Shishkevish  (1979) .  The  data  were 
obtained  from  the  recordings  of  the  Pamir-Lena  River  seismic  array 
for  earthquakes  in  the  Baikal  region,  Sinkiang,  the  Gobi  desert, 
southwestern  China,  and  the  Himalayas.  Except  for  Figures  9d 
and  9e  which  contain  propagation  paths  in  the  tectonically  ac¬ 
tive  mountain-belts  of  Central  Asia,  the  short-period  Lg/P  ratios 
are  generally  greater  than  6  for  propagation  paths  in  the  stable 
region  of  central  and  eastern  USSR. 

C.  Logarithmic  Ratios  of  Amplitude/Period  (A/T)  for  Lg  to  A/T 
for  P  vs.  Distance  ~ 

Since  the  results  from  the  last  section  did  not  take  the  epi- 
central  distances  into  consideration,  in  this  section  we  plot  the 
amplitude/period  ratios  of  Lg-  to  P-  waves  vs.  distance  in  Figures 
10  and  11  to  see  if  this  approach  may  improve  the  separation  be¬ 
tween  the  explosion-  and  earthquake-  populations.  We  have  in¬ 
cluded  the  period  of  the  observed  wave  in  the  calculation  to 
the  determination  of  body-  or  surface-wave  magnitudes,  in  the 
hope  of  reducing  the  scatter  introduced  by  the  period  differences 
of  the  observed  waves. 


For  the  western  USSR,  as  shown  in  Figure  10,  the  logarithmic 
rations  of  A/T  at  epicentral  distances  less  than  10°,  although 
sparse,  are  approximately  equal  to  zero,  i.e.  (A/T)Lg  is  approxi¬ 
mately  the  same  as  (A/T)p.  For  epicentral  distances  greater  than  10° 


Figures  9a-9e.  The  logarithmic  ratios,  A_  /A  ,  vs-  epicentral 
distance  for  earthquakes  in  (a)  the  CisbayRal^region,  (b)  Sink- 
iang,  (c)  Gobi  desert,  (d)  southwestern  China,  and  (e)  the  Him¬ 
alayas  as  recorded  by  the  seismic  stations  of  the  Pamir-Lena 
River  profile.  Solid  circles  represent  data  obtained  from  the 
short-period  SKM-3  seismographs;  solid  circles,  from  unspecified 
long-period  (probably  SKD)  seismographs;  and  shaded  regions, 
mountain  belts.  The  date  of  the  earthquakes  is  also  specified 
on  the  plot. 
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Figure  10.  Logarithmic  ratios  of  A/T  for  Lg  to  A/T  for  P 
epicentral  distance  in  western  USSR. 


1.0 


0.5 


0.0 


-0.5 


-1.0 


CENTRAL  USSR 

•-EARTHQUAKE 

x-EXPLOSION 


x 

x 


_L 


X 

X 


* 


X 

X 


10  15 

DISTANCE  (degree) 


x 

XH 

x 


*x3 


* 

X 


20 


Figure  11.  Logarithmic  ratios  of  A/T  for  Lg  to  A/T  for  P  vs 
epicentral  distance  in  central  USSR. 


the  logarithmic  ratios  show  large  scatter  and,  in  general,  have 
negative  values.  In  Figure  11,  the  logarithmic  ratios  obtained 
at  the  WWSSN  stations  KBL  and  MSH  are  shown  for  25  presumed  ex¬ 
plosions  at  the  East  Kazakh  test  site  (~  49.8°n,  78.2°E)  and  five 
Gazli  earthquakes.  Since  the  presumed  explosions  occurred  in  the 
same  region,  the  epicentral  distances  tend  to  cluster  at  16.5°- 
17.0°  for  KBL  and  19.0°-19.5°  for  MSH.  Considering  the  relative 
similarity  in  the  source  function  and  the  proximity  of  the  propa¬ 
gation  paths,  the  wide  variation  in  the  logarithmic  ratios  at  a 
single  station  is  most  striking.  In  contrast  to  the  scatter  of 
the  explosion  data,  the  logarithmic  ratios  for  the  Gazli  earth¬ 
quakes  seem  to  lie  closely  together. 

Although  the  A/T  ratios  of  Lg  to  P  in  these  two  figures  show 
large  scatter,  two  general  patterns  can  be  discerned.  Firstly, 
the  logarithmic  ratios  are  on  the  average  larger  at  near  distances 
than  farther  away.  Secondly,  the  earthquake  population  apparently 
cannot  be  separated  from  the  explosions  based  on  this  method.  The 
falloff  of  the  logarithmic  ratios  with  distance  is  probably  a  re¬ 
sult  of  differences  in  the  geometrical  spreading  and  attenuation 
characters  of  P  and  Lg  waves.  [In  a  homogeneous  sphere,  the  geo¬ 
metrical  spreading  would  introduce  a  factor  of  r  ^  to  P  waves  but 
-1/2 

only  (sin  A)  to  Lg  waves,  where  r  and  A  are  the  epicentral  distances 

in  km  and  radians,  respectively.  In  a  spherically  layered  earth  the 
spreading  factor  for  P  waves  would  depend  on  the  cross  section  of 
the  ray-pencil  at  the  source  and  the  receiver,  which  depend  on  the 
elastic  parameters  at  the  source  and  the  receiver  as  well  as  the 
layers  above  the  turning  point  of  the  ray  (cf.  Aki  and  Richards, 

1980);  but  the  spreading  factor  for  Lg  waves  would  remain  the  same. 
Similarly,  the  differences  in  the  propagation  path  would  affect 
the  attenuation  of  P-  and  Lg-  amplitudes  differently.]  If  we  as¬ 
sume  that  the  effects  of  the  geometrical  spreading  and  attenuation 
can  be  combined  at  regional  distances — again,  similar  to  the  cor¬ 
rection  term  used  in  calculating  the  body-  or  surface-wave  magnitude — 
then  the  falloff  of  the  logarithmic  ratios  with  distance  implies 
a  faster  diminution  of  Lg-  amplitude  than  P-  amplitudes:  a  factor 


that  should  be  taken  into  account  quantitatively  in  future  studies. 
Also,  since  the  usage  of  amplitude/period  ratios  do  not  improve  the 
separation  between  the  earthquakes  and  explosions  over  the  ampli¬ 
tude  ratios  noticeably,  we  would  recommend  using  the  latter  which 
is  simpler  than  the  former,  in  conjunction  with  some  other  methods. 

D.  Group  Velocity  at  Amplitude  Maxima 

Since  the  comparison  between  the  amplitudes  of  short-period 
P  waves  to  those  of  Lg  waves  at  regional  distances  did  not  prove 
*  to  be  as  useful  a  discriminant  as  the  m^  -  Mg  method  at  teleseis- 

mic  distances,  we  directed  our  efforts  to  the  search  of  depth  in¬ 
formation  from  the  coda  which  arrives  after  the  direct  S.  The 
coda  generally  contains  the  largest  amplitudes  at  regional  dis¬ 
tances  and  includes  the  multiply-reflected  S  waves  as  well  as  Lg 
waves.  The  description  in  the  first  part  of  this  final  report 
indicates  that  the  waveform  of  Lg  is  probably  a  superposition  of 
higher-mode  surface  waves  (both  Love-  and  Rayleigh- type) ,  the 
.  relative  excitation  of  which  depends  on  the  focal  depth  and  the 

source  mechanism,  modified  by  the  anelastic  and  scattering  prop¬ 
erties  of  the  propagation  path.  Also,  several  studies  (e.g. 

Knopoff  et  al.,  1973;  Panza  and  Calcagnile,  1975;  etc.)  indicate 
I  that  (i)  shallow  events  tend  to  excite  the  fundamental  and  the 

lower-order  modes  more  efficiently  than  the  higher-order  modes, 
and  (ii)  for  waves  with  the  same  period,  fundamental  and  lower- 
order  modes  generally  have  lower  group  velocities  than  the  higher- 
t  order  modes.  Thus,  if  the  effects  of  the  propagation  path  are  ne¬ 

glected,  then  shallow  events  would  tend  to  contain  larger  amplitudes 
at  lower  group-velocity  window  than  deep  events.  With  this  theor¬ 
etical  possibility  in  mind,  we  measured  the  group  velocity  at  am- 
(  plitude  maxima  as  a  function  of  distance. 

For  events  in  the  eastern  U.S.,  the  group  velocities  at  am¬ 
plitude  maxima  vs.  distance  shown  in  Figure  12  were  obtained  at 
stations  of  the  NEUSSN;  Figure  13  from  stations  of  the  V7WSSN  and  the 
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Figure  12.  Group  velocities  measured  at  amplitude  maxima  vs.  distance  for  propagation  paths 
in  the  eastern  U.S.  as  recorded  by  the  NEUSSN. 
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Figure  13.  Group  velocities  measured  at  amplitude  maxima  vs.  distance  for  propagation  paths  in  the 
eastern  U.S.  as  recorded  by  the  stations  of  the  WWSSN  and  the  University  of  Minnesota  array. 


University  of  Minnesota  array  (the  Wichita  Mountain  and  Cumber¬ 
land  Plateau  Seismic  Observatories,  abbreviated  as  WMSO  and  CPSO) ; 
and  Figure  14,  from  WWSSN  and  LRSM  stations.  Similarly,  measure¬ 
ments  in  Figures  15  and  16  were  derived  from  WWSSN  stations  for 
propagation  paths  mostly  in  the  western  and  central  USSR,  re¬ 
spectively.  An  examination  of  these  figures  indicate  that  except 
for  central  USSR,  the  group  velocities  at  the  amplitude  maxima  for 
explosions  seem  to  be  less  than  those  for  earthquakes.  This  ob¬ 
servation  would  support  our  hypothesis  if  the  earthquakes  are  lo¬ 
cated  deeper  than  the  explosions.  In  general,  the  burial  depths 
for  underground  explosions  are  less  than  1  km,  whereas  the  focal 
depths  for  the  earthquakes  used  in  this  study  are  poorly  known. 
(Using  the  travel  times  of  P  waves,  the  depth  resolution  for  shal¬ 
low  teleseismic  events  is  probably  no  better  than  +  25  km.]  Tec¬ 
tonic  considerations,  however,  can  constrain  the  focal  depths  of 
earthquakes  in  the  eastern  U.S.  and  the  western  USSR  to  be  less 
than  35  km.  The  proximity  of  the  Gazli  earthquakes  to  the  Alpine- 
Himalayan  orogenic  belt  would  probably  increase  this  uncertainty 
even  further.  The  only  conclusion  which  we  can  safely  draw  from 
the  above  discussion  is  that  the  focal  depths  of  the  earthquakes 
were  probably  deeper  than  those  of  the  explosions.  A  lack  of 
more  stringent  constraints  on  their  focal  depths  prevents  us 
from  testing  our  hypothesis  quantitatively.  For  central  USSR, 
the  large  spatial  separation  between  the  East  Kazakh  test  site  and 
Gazli  and/or  the  inclusion  of  major  tectonic  boundaries  in  the 
propagation  may  also  explain  the  disparity  in  the  observed  group 
velocities. 

E.  Energy-Ratio  Method 

During  the  initial  stage  of  the  group  velocity  study  de¬ 
scribed  above,  it  was  noted  that  for  earthquakes  in  the  eastern 
U.S.  the  energy  in  the  coda  with  the  largest  amplitude,  which 
normally  arrives  after  the  direct  S-phase,  was  distributed 
roughly  equally  about  a  group  velocity  of  3.4  km/sec.  That  is. 
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Figure  14.  Group  velocities  measured  at  amplitude  maxima  vs.  distance 
for  propagation  paths  in  the  eastern  U.S.  as  recorded  by  the  WWSSN  and 
LRSM  stations. 
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Figure  15.  Group  velocities  measured  at  amplitude  maxima  vs.  distance  for  propagation  paths 
mostly  in  the  western  USSR  as  recorded  by  the  stations  of  WWSSN. 
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Figure  16.  Group  velocities  measured  at  amplitude  maxima  vs.  distance  for 
propagation  paths  mostly  in  the  central  USSR  as  recorded  by  the  stations  of 
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about  half  of  the  energy  in  the  coda  propagated  at  a  group  velocity 
greater  than  3.4  km/sec  ,  while  the  other  half,  at  a  group  velocity 
of  less  than  3.4  km/sec.  As  a  result  of  this  observation,  two 
group-velocity  windows  were  selected:  3. 4-4.0  km/sec  and  2. 8-3. 4 
km/sec  ,  to  see  if  the  energy  distribution  in  these  two  windows 
differs  between  earthquakes  and  underground  explosions.  The  moti¬ 
vation  behind  this  approach  is  similar  to  that  of  the  previous  sec¬ 
tion  on  group  velocity,  i.e.  shallow  events  presumably  contain  more 
energy  in  the  2. 8-3. 4  km/sec  window,  during  which  the  fundamental 
and  lower-order  modes  arrive,  than  comparably-sized  deep  events. 
Thus,  instead  of  measuring  the  group  velocity  of  the  amplitude 
maxima  at  a  single  point,  the  energy-ratio  approach  averages  the 
amplitude  spectra  of  two  band-limited,  group- velocity  windows  and 
compares  them.  Moreover,  since  our  measurements  were  taken  from 
the  vertical-component  seismograms,  only  waves  of  the  Rayleigh-,  P-, 
and  SV- types  are  of  interest  to  us. 

To  quantify  the  amount  of  energy  within  each  group-velocity 
window,  we  measured  the  area  enclosed  by  the  envelope  of  the  wave¬ 
form  in  the  selected  group-velocity  windows  with  a  planimeter.  This 
technique  is  analogous  to  the  AR-method  used  by  Brune  et  al.  (1963) 
on  long-period  surface  waves.  Since  the  area  measured  is  propor¬ 
tional  to  the  energy  contained  in  the  group-velocity  window,  we 
have  designated  the  areas  in  the  3. 4-4.0  km/sec  and  2. 8-3. 4  km/sec 
windows  by  EHIGH  and  ELQW,  respectively.  The  subscripts  high  and 
low  refer  to  the  relative  group  velocity  in  the  two  windows. 

Figure  17  shows  EHIGfJ  vs.  ELQW  for  events  in  the  eastern  U.S. 
Results  from  this  figure  seem  to  indicate  that  the  underground 
nuclear  explosion,  SALMON,  contained  relatively  more  energy  in  the 
low  group- velocity  window  than  the  earthquakes.  Although  this  ob¬ 
servation  may  not  be  independent  of  the  lower  group  velocity  ob¬ 
served  for  SALMON,  the  technique  had  nonetheless  improved  the 
separation  between  the  earthquake-  and  explosion-populations  in 
this  case. 


ELOw  (2.8  -  3.4  km/sec) 


Figure  17.  Energy  in  the  3. 4-4.0  km/sec  window  (EHIGH)  vs.  energy 
in  the  2. 8- 3. 4  km/sec  window  (ET-„)  for  events  in  the  eastern  U.S. 
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Figures  18-20  show  the  ratios  of  EHIGH  to  ELQW  as  a  function 
of  epicentral  distance  for  the  eastern  U.S.,  and  the  western  and 
central  portions  of  the  USSR,  respectively.  An  inspection  of 
these  figures  show  that  the  energy  ratios  exhibit  a  clear  separ¬ 
ation  between  earthquakes  and  explosions  in  the  eastern  U.S.,  but 
not  in  the  western  USSR.  (The  plot  for  central  USSR  did  not  con¬ 
tain  any  earthquake  data;  consequently,  no  comparison  was  possible.) 
The  difference  in  the  discrimination  ability  may  be  explained  in 
several  ways.  Firstly,  the  data  from  SALMON  was  anomalous  because 
of  the  effects  of  the  unusual  burial  medium,  salt,  and  the  propa¬ 
gation  through  the  thick  sedimentary  wedge  of  the  Mississippi  Em- 
bayment.  Secondly,  the  peaceful  nuclear  explosions  (PNE’s)  con¬ 
ducted  in  western  USSR  may,  for  engineering  purposes,  have  been 
designed  or  deployed  differently  from  the  non-PNE's.  Thirdly,  the 
great-circle  paths  from  these  events  in  the  western  USSR  to  the  re¬ 
cording  stations  generally  include  one  or  several  large-scale, 
lateral  heterogeneities  (e.g.  the  Gulf  of  Finland,  Gulf  of  Bothnia, 
Baltic  Sea,  or  the  Alpine-Himalayan  belt)  which  may  affect  the 
energy  distribution  in  the  coda  by  frequency-dependent  absorption, 
scattering,  and  changes  in  group  velocity.  Lastly,  the  selected 
group-velocity  windows  may  have  to  be  modified  in  different  regions 
to  optimize  the  potential  of  extracting  depth  information  from  the 
Lg  coda.  Other  or  a  combination  of  these  explanations  is,  of  course, 
also  quite  possible. 


This  paragraph  expands  on  the  optimization  of  the  energy-ratio 
method  mentioned  above.  According  to  the  study  of  Herrmann  (1974) , 
the  excitation  function  of  higher-mode  surface  waves  depends  primar¬ 
ily  on  the  focal  depth.  Thus,  if  we  know  the  average  structure  of 
the  source  and  the  propagation  path,  then  we  can  calculate  the  rela¬ 
tive  importance  of  the  lower-order  to  higher-order  modes  for  sources 
at  a  certain  depth,  say,  5  km.  Since  we  can  also  calculate  the  group- 
velocity  curves  for  the  higher-mode  surface  waves,  we  should  be  able 
to  define  a  threshold  group  velocity  at  which  the  difference  in  the 
energy  ratio  for  the  given  hypocentral  depth  is  maximized.  This  ap¬ 
proach  will  be  pursued  in  future  studies  on  regional-distance  dis¬ 
crimination. 


I 


LU 

<  z 

Z>  o 
o 

x  o 
K-  -J 
or  cl 

<X  X  A 
LU  LU  U 

I  ! 

O  X 


o  o 


*o 

o 


X 

X  X 


©»>o 


eo  x 


XOKX 


w 

o 

• 

10 

m 

0 

• 

-H 

D 

+J 

(0 

C 

« 

u 

0) 

4J 

CO 

to 

DISTANCE  (degree) 


T 


T 


1.5  h  x— EXPLOSION 


1.0  b 


0.5 


0.0 


CENTRAL  USSR 


X 

* 

X 


10 


X 

X 


$ 


15  20 

DISTANCE  (degree) 


Figure  20.  Ratios  of  to  ELQW  as  a  function  of  epicentral 


distance  for  propagation  paths  in  the  central  USSR. 
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The  various  discrimination  methods  discussed  in  this  section 
all  seem  to  indicate  that  significant  discrepancies  exist  in  the 
results  derived  from  the  different  regions.  Differences  in  the 
propagation  characteristics  of  the  various  regions  can  explain 
most  of  the  observed  discrepancies.  Thus,  we  will  need  to  improve 
our  understanding  on  the  propagation  characteristics  of  seismic 
waves  on  a  regional  basis  so  that  we  may  (i)  assess  the  feasibility 
of  the  methods  discussed  above  with  more  confidence,  and  (ii)  devise 
some  other  discrimination  methods. 


Part  III.  Preliminary  Studies 


As  part  of  our  program  on  the  seismic  wave  propagation  at  re¬ 
gional  distances  we  have  initiated  several  studies,  the  prelimin¬ 
ary  results  of  which  are  reported  below. 

A.  Attenuation  of  Lg  Waves 

Studies  on  the  attenuation  of  Lg  waves  were  carried  out  for 
the  eastern  U.S.,  as  well  as  western  and  central  USSR.  In  the 
eastern  U.S.,  readings  on  the  amplitude  and  the  period  of  Lg  were 
made  from  the  short-period  seismograms  of  WWSSN  and  NEUSSN.  The 
amplitudes  were  normalized  relative  to  the  assigned  magnitude  (by 
USGS,  NOAA,  or  St.  Louis  University)  of  the  event  and  the  ratios 
of  normalized  amplitude  to  period  vs.  epicentral  distances  were  then 
plotted  on  Figure  21.  Since  the  assiqned  magnitude  was  probably 
derived  from  averaging  a  limited  number  of  readings,  measurements 
at  different  stations  would  inevitably  deviate  from  this  mean. 
Moreover,  since  these  deviations  are  propagated  into  the  normal¬ 
ization  process,  the  scaled  amplitudes  would  not  only  be  subject 
to  the  modulating  effects  at  the  recording  site  but  also  to  those 
from  which  the  assigned  magnitude  was  based  on.  If  our  normal- 
ized-amplitude/period  ratios  can  approximate  statistically  the  un¬ 
biased  values,  then  the  data  shown  in  Figure  21  suggests  a  slightly 
higher  attenuation  rate  for  the  eastern  U.S.  than  that  derived  by 
Nuttli  (1973)  for  the  central  U.S. 

The  data  from  the  western  and  central  USSR  were  obtained  some¬ 
what  differently  since  most  of  the  events  were  presumed  explosions. 
The  procedure  used  to  obtain  the  amplitude  vs.  epicentral  distance 
relation  is  as  follows:  (i)  The  estimated  yields  from  Dahlman  and 
Israelson  (1977)  were  converted  to  body-wave  magnitudes  via  the  re¬ 
lation 

m^  =  0.93  iog^Q  Y  +  3.49. 

This  empirical  relation  was  derived  by  Ericsson  (1971)  for  data 
from  the  NTS  explosions.  Ericsson  claimed  that  this  relation  is 


A/T  (/j./ sec) 


Figure  21.  Ratios  of  normalized  amplitude  to  period  of  Lg  waves 
vs.  epicentral  distance  in  the  eastern  U.S.  See  text  for  a  dis¬ 
cussion  of  the  normalization  procedure. 
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also  representative  for  tests  in  the  USSR.  The  validity  of  this 
claim,  however,  remains  to  be  demonstrated.  (ii)  We  then  normal¬ 
ized  the  observed  amplitudes  corresponding  to  an  m^  =  4.4  event 
by  using  the  following  formula 

Normalized  Amplitude  =  Observed  Amplitude  X  10 
(This  second  step  is  similar  to  the  procedure  used  to  normalize 
the  events  in  eastern  U.S.).  Since  the  magnitude  used  in  the 
above  calculation  was  estimated,  the  uncertainties  in  the  normal¬ 
ized  amplitudes  would  probably  exceed  the  ones  for  earthquakes. 

The  procedure  nevertheless  provides  a  first-order  estimate  for  the 
attenuation  of  Lg  waves. 

Explosion  data  from  the  western  and  central  USSR,  earthquake 
data  from  the  western  USSR,  as  well  as  the  normalized  amplitude- 
mean  at  500,  1000,  1500,  and  2000  km  taken  from  Antonova  et  al. 

(1978)  in  Figure  22.  The  earthquake  data  from  the  western  USSR, 
was  normalized  similar  to  that  from  the  eastern  U.S.  The  re¬ 
sults  from  Antonova  et  al.  (1978)  were  recorded  at  the  Pamir-Lena 
River  seismic  array  for  earthquakes  in  the  Central  Asia.  Antonova 
et  al.  concluded  that  (i)  at  epicentral  distances  less  than  700  km, 
the  amplitudes  of  Lg  are  proportional  to  ‘  ,  and  (ii)  the  ex¬ 

ponent  decreases,  i.e.  becomes  more  negative,  as  distance  increases 
such  that  at  2000  km  the  exponent  is  approximately  between  -2.2  and  -2.5. 
Despite  the  large  scatter,  our  data  is  not  inconsistent  with 
curves  having  slopes  between  -2  and  -3  at  these  epicentral  dis¬ 
tances.  In  comparison  to  the  results  from  the  eastern  U.S.,  the 
attenuation  rate  in  the  western  and  central  USSR  appears  rather 
high.  But,  if  we  take  the  propagation  paths,  which  straddle  one 
or  several  major  tectonic  boundaries  (most  of  the  Soviet  data  was 
derived  from  earthquakes  outside  the  USSR,  while  our  data  was  ob¬ 
tained  from  presumed  Soviet  explosions  as  recorded  outside  the 
USSR) ,  into  consideration  then  the  higher  attenuation  rates  are 
not  unreasonable. 

B.  Attenuation  and  Magnitude-Scale  for  Intermediate-Period 
Rayleigh  Waves  ~ 


For  propagation  paths  in  the  eastern  North  America,  Rayleigh 


waves  with  periods  between  8  to  14  seconds  are  the  most  prominent 
feature  on  the  long-period  seismograms  of  WWSSN  (To=15  sec,  Tg=100  sec) 

(a)  Attenuation 


To  measure  the  anelastic  properties  of  intermediate-period 
Rayleigh  waves,  we  measured  the  amplitudes  and  periods  at  the 
amplitude  maxima  and  plotted  the  ratios  of  amplitude  to  period 
as  a  function  of  epicentral  distances  in  Figure  23.  In  a  study 
on  the  surface-wave  attenuation  of  central  U.S.,  Nuttli  (1973) 
showed  that  in  the  distance  range  2°  to  20°  the  falloff  of  wave 
amplitude  with  propagation  distance,  due  to  the  effects  of  geo¬ 
metrical  spreading  and  anelastic  attenuation,  can  be  approxi¬ 
mated  by  a  straight  line  on  a  log-log  plot.  Following  Nuttli' s 
example,  we  found  that  the  data  in  Figure  23  can  be  fitted  by  a 
straight  line  with  a  slope  of  -1.66,  which  corresponds  to  an  at¬ 
tenuation  coefficient  of  0.10  deg  1.  This  attenuation  rate  is 
the  same  as  that  derived  by  Nuttli  for  the  central  U.S.  but  dif¬ 
ferent  from  those  of  Basham  (1971)  and  Evernden  et  al.  (1971)  . 

We  concur  with  Nuttli 's  (1973)  interpretation  that  the  discrep¬ 
ancy  arises  from  the  phase  Rg,  instead  of  the  fundamental-mode 
Rayleigh  waves,  measured  by  Basham  and  Evernden  et  al.  (Nuttli 
also  notes  that  Rg  is  "...prominent  on  the  seismograms  of  North 
American  stations  for  earthquakes  or  underground  explosions  in 
the  western  United  States.  However,  it  is  not  well  developed 
for  earthquakes  in  the  central  United  States  recorded  at  eastern 
stations...").  The  magnitude- scale  derived  by  Vanek  et  al.  (1962) 
indicates  that  the  amplitudes  of  surface  waves,  at  periods  near 
10  sec  ,  are  attenuated  at  the  same  rate  as  in  the  eastern  and 
central  U.S.  This  observation  together  with  our  findings  in  Part  II 
on  Lg-propagation  in  the  U.S.  and  the  USSR  seem  to  imply  a  close 
similarity  between  the  crustal  structure  in  the  eastern  U.S.  and 
that  in  the  eastern  and  central  USSR. 

(b)  Magnitude  Scale 

Based  on  the  attenuation  rate  derived  above,  we  obtained  a 
magnitude-scale  formula  for  regional  distances 


Figure  23.  Ratios  of  amplitude  to  period  of  intermediate-period 
(8-14  sec)  Rayleigh  waves  vs.  epicentral  distance  in  the  eastern 
U.S.  The  dotted  line  represents  an  approximation  of  7T=0.1  deg 
in  the  distance  range  of  2°  to  20°. 
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MS  =  log10  (Amax/T)  +  1 -66  log10A  +  2,60 
This  formula  was  slightly  modified  from  Equation  (4)  of  Nuttli 
(1973) ,  where  he  used  (A/T)  and  wave  periods  of  3-12  seconds 

IHclX 

instead  of  (A  /T)  and  8-14  seconds  used  in  this  study.  In  both 
max 

cases,  the  range  of  applicability  is  between  2°  and  20°  in  eastern 
North  America. 

C.  Intermediate-Period  M„  vs.  M„  in  Eastern  and  Central  U.S. 

- - - s - Lg - 

Having  determined  the  magnitude-scale  formulae  for  Lg  waves 
(MLg  at  0. 3-1.0  sec  )  and  intermediate-period  Rayleigh  waves  (Mg 
at  8-14  sec  )  appropriate  for  the  eastern  and  central  U.S.,  we 
became  interested  in  investigating  (i)  the  relationship  between 
Mg  and  M^,  and  (ii)  the  possibility  of  using  them  as  a  discrim¬ 
inant.  The  magnitudes  Mg  and  ML  for  four  earthquakes  and  one 
underground  nuclear  explosion  (SALMON)  in  the  eastern  U.S.  were 
measured  and  plotted  in  Figure  24;  also  shown  in  this  figure  are 
the  data  points  for  four  central  U.S.  earthquakes  taken  from  a 
study  by  Nuttli  (1973).  An  inspection  of  this  figure  shows  that 
(i)  all  the  data  points  can  be  approximated  by  a  linear  relation 
of  the  form:  Mg  =  1.69  M^  -  4.08;  and  (ii)  the  explosion  data 
cannot  be  discriminated  from  the  earthquake  population.  Since 
there  is  only  one  explosion  used  in  this  preliminary  study,  we 
would  like  to  examine  more  data  from  eastern  and  central  U.S. 
explosions  (e.g.  RULISON,  GAS  BUGGY,  RIO  BLANCO,  etc.)  in  the  future. 

D.  Usefulness  of  High-Frequency  Waves  at  Regional  Distances 

Although  many  of  the  NEUSSN  stations  operate  with  peak  mag¬ 
nifications  in  the  10-20  Hz  range,  the  seismograms  for  earthquakes 
within  regional  distances,  as  examined  by  us,  did  not  show  fre¬ 
quencies  higher  than  5  Hz.  This  observation  is  quite  different 
from  the  efficient  propagation  of  Lg  and  intermediate-period  sur¬ 
face  waves  in  the  eastern  and  central  U.S.;  consequently,  we  sus¬ 
pect  that  waves  at  frequencies  higher  than  5  Hz  are  attenuated 
rapidly  by  scattering  at  small-scale  heterogeneities  so  that  they 
may  not  be  very  useful  at  regional  distances  in  certain  regions. 


EASTERN 

U.S.  o 

/  Ms  =  1.69Ml§ 

/ 

o  /  •  0 

/ 

/ 

/ 

/ 

S* 

/ 

/ 

s 

•  This  study 

/ 

o  Nuttli  (1973) 

A  • 

x  SALMON 

Figure  24.  Mg  (from  8-13  sec  Rayleigh  waves)  vs.  M  (from  0.3 
1.0  sec  Lg  waves)  for  events  in  the  eastern  and  central  U.S. 


Part  IV. 


Magnitude-Yield  Relation  and  Others 


An  accurate  determination  of  the  magnitude-yield  relation  is 
an  important  geophysical  problem.  Aside  from  its  obvious  appli¬ 
cation  for  estimating  the  yield  of  unknown  nuclear  tests  by  meas¬ 
uring  the  amplitudes  of  the  observed  seismic  waves,  a  well-deter¬ 
mined  magnitude-yield  relation  may  become  one  of  the  most  useful 
tools  for  calibrating  the  seismic  energy  (especially  at  short 
periods)  radiated  by  earthquakes.  The  task  of  casting  this  rela¬ 
tion  into  a  well  defined  form,  however,  is  not  an  easy  one.  Dif¬ 
ficulties  can  be  traced  to  both  the  magnitude  and  the  yield  ends 
of  the  relation.  Below  we  will  describe  some  of  the  difficulties 
involved. 

The  amplitudes  of  the  observed  seismic  waves  can  be  signifi¬ 
cantly  affected  by  several  factors,  such  as  (i)  the  medium  and 
the  burial  depth  of  the  source,  (ii)  the  degree  of  seismic  coup¬ 
ling  between  the  source  and  the  surrounding  medium,  and  (iii)  the 
local  structures  beneath  the  source  and  the  receivers.  The  first 
and  third  factors  have  plagued  seismologists  for  years,  but  these 
problems  are  currently  being  solved.  To  our  knowledge,  the  second 
factor  has  not  been  studied  extensively,  its  effects  are  therefore 
not  well  understood. 

Several  investigators  have  attempted  to  establish  the  magni¬ 
tude-yield  relation  b&sed  on  magnitudes  that  are  determined  from 
local/regional  networks  and/or  a  relatively  small  number  of  events. 

In  view  of  the  lack  of  completeness  of  these  studies  and  the  im¬ 
portance  of  this  problem,  we  have  decided  to  (i)  undertake  a  com¬ 
prehensive  compilation  of  available  published  results  that  are 
relevant  to  the  problem  of  yield-estimation,  (ii)  present  the  re¬ 
sults  from  our  compilation  in  a  useful  form,  and  (iii)  improve  the 
determination  of  body-wave  magnitudes,  in  a  statistical  sense,  by 
increasing  the  number  of  amplitude  measurements  at  various  epi- 
central  distances.  [ISC  determines  its  body-wave  magnitudes  only 
if  three  or  more  stations  report  their  amplitudes.  It  then  applies  the 


*  * 


unified  magnitude  of  Gutenberg  (1956)  to  the  amplitudes  to  determine 
the  m^.  Few  stations,  however,  have  the  habit  of  reporting  their 
amplitudes  to  the  ISC.] 

Data 

Because  of  the  large  number  ( *>  400)  of  nuclear  tests  in  the 
United  States  and  the  Soviet  Union,  we  have  limited  most  of  our 
data  base  to  those  underground  nuclear  explosions  for  which  re¬ 
ports  on  their  estimated  yield  exist.  The  U.S.  data  used  is  de¬ 
rived  from  Springer  and  Kinnaman  (1971,  1975) ,  and  the  Soviet  data, 
from  Bolt  (1976)  and  Dahlman  and  Israelson  (1977) .  The  magnitude 
determinations  used  are  from  Bolt  (1976)  and  the  International 
Seismological  Centre  (ISC)  Bulletins.  There  are  some  doubts  con¬ 
cerning  the  source  reference  of  the  estimated  yield  for  the  Soviet 
tests,  compiled  by  Dahlman  and  Israelson,  as  well  as  the  magni¬ 
tude  of  the  Soviet  tests  as  reported  by  Bolt;  we  are  in  the# pro¬ 
cess  of  uncovering  these  uncertainties. 

Table  VI  represents  a  compilation  of  the  U.S.  explosion  data 
used  in  this  report.  The  table  contains  the  name,  data,  origin 
time,  location,  and  burial  depth  of  the  event;  it  also  describes 
the  rock-type  surrounding  the  buried  source  (e.g.  tuff,  alluvium, 
rhyolite,  etc.),  the  dimensions  (volume,  diameter,  and  height)  of 
the  collapse  cavity,  the  body-wave  magnitude  (ISC) ,  and  the  an¬ 
nounced  or  estimated  yield.  Except  for  the  magnitude,  all  the 
information  was  provided  to  Springer  and  Kinnaman  (1971,  1975) 
by  the  U.S.  Atomic  Energy  Commission  (AEC) .  A  compilation  of  the 
available  Soviet  data  is  outlined  in  Table  VII.  A  compilation  of  the 
of  the  date,  computed  origin  time  and  location  (Bolt,  1976) ,  the 
body-wave  magnitudes  (from  ISC  and  Bolt's  compilation),  and  the 
estimated  yield  for  these  events  (Dahlman  and  Israelson,  1977) . 

Based  on  the  compilations  in  Table  VI  and  VII,  we  have  made  the 
following  plots: 

From  the  Soviet  data:  (ISC  and  Bolt's)  vs.  estimated 

yield  (Figure  25  and  26,  respectively) 


Figure  25.  Body-wave  magnitude  (ISC)  vs.  yield  for  events  in  the 
USSR.  The  dashed  line,  m.  =  0.75  log.gY  +  4.345,  represents  our 

inqDrelation  betw« 


preliminary,  best-fitting 


Eween  these  two  quantities. 


From  the  U.S.  data: 


a.  m^  (ISC)  vs.  estimated  yield 
(Figure  27) 

b.  volume  of  collapse  vs.  esti¬ 
mated  yield  (Figure  28) 

c.  diameter  and  height  of  collapse 
center  vs.  estimated  yield 
(Figures  29  and  30  respectively) 

d.  volume  of  collapse  vs.  depth  of 
burial  (Figure  31) 

Information  on  the  locality  and  the  rock-type  of  the  test-site 
are  also  included  whenever  available. 

Results  and  Discussion 


A  comparison  between  the  empirically  determined  and  computed 
magnitude-yield  relations  in  different  media  (cf.  Figures  7-8  of 
Bolt,  1976)  and  the  data  points  in  Figure  25  and  27  shows  that 
the  U.S.  data  can  be  approximated  closely  by  the  curve  for  granite, 
whereas  the  Soviet  data  lies  roughly  between  the  curves  for  granite 
and  water.  Body-wave  magnitudes  taken  from  Bolt,  on  the  other  hand , 
show  larger  scatter  than  m^  (ISC)  when  plotted  as  a  function  of 
estimated  yield  (Figures  25  and  26) .  There  is  some  indication  that 

(i)  events  in  the  East  Kazakhstan  are  more  efficient  in  generating 
seismic  waves  than  the  other  test  sites  of  the  Soviet  Union,  and 

(ii)  events  situated  in  tuff  and  rhyolite  generate  waves  more 
efficiently  than  those  located  in  alluvium  at  the  Nevada  Test  Site 
(NTS) . 

In  plotting  the  collapse  volume  vs.  the  estimated  yield 
(Figure  28) ,  we  divided  the  data  into  three  groups:  the  first 
two  groups  (open  and  closed  symbols)  refer  to  events  presented 
in  Figure  27,  while  the  third  group  (semi-filled  symbols)  con¬ 
sists  of  events  that  contain  information  on  the  collapse  volume 
and  the  estimated  yield  but  not  on  the  body-wave  magnitude.  The 
first  two  groups  are  divided,  somewhat  arbitrarily,  into  normal 


Volume  (yd3) 


Figure  28.  Volume  of  collapsed  crater  vs.  estimated  yield  for  events 
in  the  U.S.  Open  and  filled  symbols  refer  to  events  presented  in 
Figure  3,  whereas  semi-filled  symbols  refer  to  events  that  contain  in¬ 
formation  on  the  collapse  volume  and  the  estimated  yield  but  not  on 
the  body-wave  magnitude  and  therefore  not  plotted  in  Figure  3.  Filled 
symbols  denote  normal  (N)  events  which  lie  closely  together  as  a  group; 
open  symbols,  anomalous  (A)  events  which  appear  to  have  unusually  small 
collapse  volumes  for  their  estimated  yields. 


Y  (  kt ) 


Figure  29.  Diameter  of  the  collapsed  crater  vs.  estimated  yield  for 
events  in  the  U.S.  Except  for  the  semi-filled  symbols,  the  legends  are 
similar  to  those  in  Figure  4 . 


c 


Y(kt) 


H  ( ft ): 


Figure  30.  Height  of  the  collapsed  crater  vs.  estimated  yield  for 
events  in  the  U.S.  Except  for  the  semi-filled  symbols,  the  legends 
are  similar  to  those  in  Figure  4 . 


75 


(closed  symbols)  and  anomalous  (open  symbols)  events.  The  normal 
events  lie  closely  together  as  a  group,  while  the  anomalous  events 
appear  to  have  unusually  small  collapse  volumes  for  their  esti¬ 
mated  yields.  Figures  29  and  30  (the  diameter  and  depth,  respec¬ 
tively,  of  the  collapse  crater  vs.  estimated  yield)  were  plotted 
from  the  same  data  set.  It  is  quite  interesting  that  except  for 
the  anomalous  events,  the  diameter  of  the  collapse  crater  can  be 
approximated  as  being  linearly  proportional  to  the  logarithm  of 
the  yield;  the  height  of  the  crater,  however,  appears  to  be  inde¬ 
pendent  of  the  yield.  Figure  31,  which  relates  the  collapse 
volume  to  the  burial  depth,  is  composed  of  the  events  found  in 
Figure  29  (or  30  as  well  as  events  without  reports  on  their  mag¬ 
nitude  and  yield.  This  figure  seems  to  indicate  three  depth-de¬ 
pendent  distributions:  (i)  the  volume  of  collapse  is  independent 
of  burial  depth  when  the  latter  is  less  than  about  900  feet,  (ii) 
at  depths  between  900  and  2500  feet,  the  logarithm  of  the  collapse 
volume  is  approximately  linearly  proportional  to  the  burial  depth, 
and  (iii)  for  the  three  events  at  deeper  than  4000  feet,  the  vol¬ 
ume  of  collapse  is  again  unpredictable.  A  cautionary  remark  is 
deemed  necessary  at  this  point:  the  burial  depth  of  the  test 
charge  is  usually  commensurate  with  its  size;  consequently,  the 
collapse  volume  is  probably  a  complex  function  of  the  local  rock 
type,  burial  depth,  and  the  actual  yield. 
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Device  Epicenter  Collapse  Crater 

Fear  Oat*  Shot  Time  Depth  (ft)  Latitude  (°N)  lonptitude  (°W)  Medium  Volume  (yd*)  01  an  *  Ht.  (ft)  Announced  Estimated  (ISC)  Type  Name 


1961 

1201 

2X500 

1191 

37.05 

116.03 

A 

1962 

0109 

161000 

992 

37.06 

116.04 

A 

0118 

180000 

856 

37.05 

116.03 

A 

0110 

180000 

1191 

37.05 

116.04 

A 

0208 

180000 

595 

37.13 

116.05 

A 

0209 

163000 

786 

37.04 

116.04 

A 

0219 

163000 

492 

37.05 

116.03 

A 

0219 

175000 

696 

37.13 

116.04 

T 

0221 

180000 

1000 

37.13 

116.05 

A 

0101 

191000 

1191 

37.04 

116.03 

A 

0M5 

181500 

110 

37.11 

116.37 

Basalt 

0108 

180000 

841 

37.12 

116.05 

A 

011S 

163000 

784 

37.04 

116.03 

A 

0128 

180000 

614 

37.12 

116.03 

7 

0311 

180000 

448 

37.05 

116.04 

A 

0«  OS 

180000 

856 

37.04 

116.02 

A 

0406 

180000 

766 

37.12 

116.04 

A 

0421 

184000 

634 

37.12 

116.03 

A 

0427 

180000 

714 

37.12 

116.04 

A 

0S07 

193300 

848 

37.05 

116.03 

A 

0512 

190000 

1424 

37.07 

116.03 

T 

0519 

150000 

714 

37.12 

116.05 

A 

0525 

150000 

632 

37.13 

116.05 

A 

0601 

170000 

539 

37.05 

116.03 

A 

0606 

170000 

860 

37.05 

116.04 

A 

0621 

170000 

854 

37.04 

116.03 

A 

0627 

180000 

1340 

37.04 

116.04 

A 

0630 

213000 

489 

37.12 

116.05 

A 

0706 

170000 

635 

37.18 

116.05 

A 

0713 

160000 

1356 

37.06 

116.03 

A 

0727 

210000 

491 

37.13 

116.06 

A 

0824 

150000 

744 

37.12 

116.04 

A 

0824 

170000 

676 

37.05 

116.02 

A 

0914 

171000 

711 

37.04 

116.02 

A 

0920 

170000 

792 

37.06 

116.03 

A 

1005 

170000 

1622 

37.14 

116.05 

T 

1019 

180000 

792 

17.04 

116.02 

A 

1027 

150000 

1048 

37.15 

116.05 

A 

1127 

180000 

747 

37.12 

116.03 

T 

1207 

190000 

993 

37.05 

116.03 

A 

1212 

184500 

761 

37.05 

116.02 

A 

1961 

0208 

160000 

994 

37.15 

116.05 

A 

0208 

183000 

856 

37.05 

116.02 

A 

0129 

154900 

917 

37.04 

116.02 

A 

0405 

175200 

791 

37.04 

116.02 

A 

0522 

154000 

1289 

37.11 

116.04 

T 

0614 

141000 

642 

37.05 

116.02 

A 

0812 

234500 

992 

37.04 

116.02 

A 

0815 

1X000 

738 

37.15 

116.08 

A 

1011 

140000 

357 

37.04 

116.02 

A 

1114 

160000 

854 

37.04 

116.02 

A 

1122 

173000 

987 

37.12 

116.04 

A 

1204 

163830 

860 

37.04 

116.03 

A 

1212 

160200 

540 

37.13 

116.04 

A 

1964 

0116 

160000 

1610 

37.14 

116.05 

7 

0123 

160000 

868 

37.13 

116.04 

7 

0220 

153000 

1616 

37.15 

116.04 

7 

0111 

160200 

376 

37.05 

116.01 

A 

0414 

144000 

668 

37.13 

116.03 

7 

0415 

143000 

491 

37.04 

116.02 

A 

0424 

201000 

1663 

37.15 

116.05 

A 

0429 

204700 

859 

37.04 

116.03 

A 

0514 

144000 

536 

37.12 

116.04 

A 

0515 

161500 

792 

37.04 

116.01 

A 

2.16  COS 

584  X 

50 

13 

0 

Fisher 

9  E03 

380  X 

7 

4.7 

0 

Stoat 

1  EOS 

452  X 

42 

5.8 

0 

Apouti 

1.6  E04 

570  X 

14 

0 

Oo  mouse 

7.35  E04 

406  X 

40 

3.1 

0 

Stillwater 

8  £04 

446  X 

37 

6.5 

0 

Armadillo 

3  E04 

314  X 

42 

1.8 

0 

Chinchilla 

9.54  E03 

296  X 

11 

0 

Codsaw 

7.39  E04 

464  X 

36 

12 

0 

Cimarron 

1.8  E05 

620  X 

40 

0 

Pampas 

265  X 

at 

0.43 

Danny  Boy 

4.31  E04 

350  X 

31 

8.4 

0 

Bratos 

1.6  E05 

484  X 

67 

0 

Hopnose 

2.65  E04 

294  X 

20 

3.4 

0 

HOOSiC 

1.7  E04 

260  X 

24 

0 

Chinchilla 

3  E05 

520  X 

97 

11 

0 

Dormouse 

Prime 

2.47  EOS 

500  X 

70 

0 

Pasaic 

1.78  E04 

310  X 

17 

0 

Dead 

1.09  £05 

394  X 

72 

0 

Black 

1.1  EOS 

454  X 

62 

0 

Pace 

4  E05 

820  X 

75 

36 

0 

Aardvark 

9.8  £03 

250  X 

13 

0 

Eel 

1.06  £05 

460  X 

51 

0 

White 

1.6  E04 

300  X 

26 

0 

Raccoon 

9.7  £04 

530  X 

44 

0 

Pack  rat 

2.6  COS 

S58  X 

92 

0 

Daman  1 

7.5  EOS 

896  X 

103 

46 

0 

Haymaker 

8.7  E04 

360  X 

70 

0 

Sacramento 

6.5  £06 

1280  X 

320 

100 

0 

Sedan 

2.24  EOS 

680  X 

50 

0 

Merrimac 

6.5  E03 

340  X 

58 

0 

Wichita 

2.26  £05 

500  X 

79 

0 

fork 

7  £04 

400  X 

44 

0 

BoDac 

2.4  E05 

510  X 

97 

0 

Hyraa 

3  E05 

560  I 

92 

0 

Peba 

8.76  EOS 

850  X 

125 

110 

0 

Missis¬ 

sippi 

3.9  EOS 

606  X 

125 

- 

0 

Bandicoot 

3.9  £04 

400  X 

21 

0 

Santee 

1.25  E05 

460  X 

61 

0 

Anacostla 

2.5  E05 

510  X 

81 

0 

lendrac 

1.4  EOS 

560  X 

50 

0 

NunOat 

1.52  E05 

446  X 

73 

0 

Casselman 

2.1  £05 

514  X 

78 

0 

Acushi 

1.7  E05 

550  X 

48 

0 

Gerbil 

2.1  £05 

474  X 

114 

0 

Ferret 

Prime 

7.35  £05 

852  X 

88 

0 

Slones 

6  £04 

302  X 

58 

0 

Mataco 

1.4  E05 

496  X 

61 

.  c 

Pekan 

4.0  £04 

300  X 

40 

0 

Satsop 

2.0  £05 

480  X 

100 

0 

Grunt  on 

2.1  COS 

520  X 

69 

0 

Anchovy 

6.4  £04 

400  X 

42 

0 

Greys 

1.9  £05 

490  X 

78 

0 

Sardine 

9.3  E04 

400  X 

60 

0 

Cafle 

1.59  £06 

1040  X 

125 

19 

5.2 

N 

Fore 

1.01  C05 

960  X 

39 

0 

Oconto 

1.03  £06 

920  X 

126 

24 

5.1 

H 

Klickitat 

2.9  £04 

240  X 

46 

0 

Pike 

1.89  £04 

295  X 

22 

0 

Hook 

1.4  £04 

236  X 

19 

0 

Stwr9*on 

1.44  £06 

1265  X 

93 

100 

5.2 

0 

Turf 

1.8  £05 

480  X 

72 

15 

4.1 

N 

Pipe  Fish 

8.5  £04 

400  X 

55 

0 

Sackiwtnf 

6.3  E04 

416  X 

28 

0 

Minnow 

77 


Table  VI 

U.S.  Underground  Nucletr  Caplostsns 


Tt*r 

Oat* 

Shot  Tin* 

Device 

Oepth  (ft) 

Epicenter 

Latitude  (°N)  longtltude  (°WJ 

Medfui 

col  lapse  Crater 

Mum  (yd3)  Olaa  2  Ht.  (ft) 

Announced 

Estlaeted 

Type 

ft 

1964 

0625 

133000 

673 

37.11 

116.03 

T 

7.93  204 

416  X  47 

• 

0 

Fjd c 

coot. 

0630 

133300 

847 

37.17 

116.06 

A 

7.88  E04 

347  X  90 

9 

A 

Out) 

0716 

131500 

1277 

37.18 

116.04 

T 

1.79  205 

536  X  65 

0 

M 

0904 

181500 

856 

37.02 

116.02 

A 

1.6  205 

450  X  74 

12 

ft 

Go* nay 

1002 

200300 

1484 

37.08 

116.01 

T 

12 

4.0 

N 

Auk 

1009 

140000 

1325 

37.15 

116.08 

A 

3.89  204 

475  I  72 

38 

30 

4.0 

A 

Par 

1016 

155930 

849 

37.04 

116.02 

A 

1.2  205 

472  X  50 

0 

Barbel 

1105 

150000 

1319 

37.17 

116.07 

Ooloaf te 

12 

9 

4.8 

UAfldwjr 

1205 

211500 

1323 

37.11 

116.05 

T 

3.2  205 

740  X  60 

10 

4.8 

N 

Crepe 

1216 

200000 

692 

37.03 

116.01 

A 

1.4  204 

260  X  18 

1.2 

n 

Pjrrot 

1216 

201000 

498 

37.18 

116.07 

T 

1.9  204 

254  I  21 

2.7 

N 

NmJ pack 

1965 

0114 

160000 

706 

37.12 

116.02 

T 

1.28  205 

450  X  60 

0 

Uuol 

0204 

153000 

762 

37.13 

116.06 

A 

5.41  204 

360  X  45 

0 

Cat  n.*r* 

0216 

173000 

972 

37.05 

116.02 

A 

1.7  £05 

510  X  55 

10 

it 

*rr)  in 

0218 

161847 

588 

36.82 

115.95 

A 

1.66.  EOS 

300  X  100 

0 

Wiifadune 

0303 

191300 

2459 

37.06 

116.04 

T 

6$ 

Wagtail 

0326 

153408 

1761 

37.15 

116.04 

T 

3.B  205 

920  I  185 

35 

N 

Cup 

0405 

210000 

1466 

37.03 

116.02 

A 

1.9  204 

450  X  8 

0 

kestrel 

0414 

131400 

280 

37.28 

116.52 

ft 

4.69  204 

338  X  79 

4.3 

0 

Palanquin 

0421 

220000 

1000 

37.01 

116. ZO 

T 

8 

S.O 

GuiMlrop 

0507 

154711 

624 

37.14 

116.07 

A 

6.25  204 

385  X  30 

0 

Tee 

0S21 

130852 

922 

37.12 

.  116.03 

T 

8.92  205 

630  X  148 

0 

Tweed 

0611 

194500 

593 

37.04 

116.02 

A 

9.5  203 

290  X  17 

1.2 

0 

Petrel 

0723 

170000 

1741 

37.10 

116.03 

T 

7.9  205 

1055  X  77 

60 

S.4 

II 

Bronze 

0306 

172330 

1053 

37.02 

116.04 

A 

2.2  205 

526  X  77 

28 

ft 

ftau** 

0901 

200300 

990 

37.02 

116.01 

A 

2.1  205 

506  X  72 

12 

4.2 

n 

Screamer 

0910 

171200 

1494 

37.08 

116.02 

T 

3.5  £05 

978  X  40 

0 

Charcoal 

1112 

180000 

791 

37.05 

116.02 

A 

8.6  204 

456  X  39 

0 

Sepia 

1203 

151302 

2236 

37.16 

116.05 

T 

2.33  206 

800  I  100 

ft 

CcOuroy 

1216 

191500 

1642 

37.07 

116.03 

r 

4.2  £05 

1284  X  44 

36 

5.3 

ft 

32 

5.2 

LjAfUiiU 

1966 

0118 

183500 

1842 

37.09 

116.02 

T 

0121 

182800 

1093 

37.03 

116.02 

A 

3.8  204 

464  X  16 

0 

Dovekie 

0203 

181737 

836 

37.13 

116.07 

A 

2.57  204 

260  X  27 

0 

Plata  11 

0224 

155507 

2204 

37.27 

116.43 

T 

16 

7 

5.0 

Re* 

0307 

164100 

642 

37.04 

116.03 

A 

8.07  204 

408  X  58 

0 

7  Inf  Out 

0318 

190000 

1092 

37.01 

116.01 

T 

8.3  204 

458  X  39 

0 

Purple 

0406 

135717 

739 

37.14 

116.14 

T 

1.25  £05 

386  X  85 

5 

4.4 

N 

Slut* 

0407 

222730 

742 

37.02 

115.99 

I 

2.5  204 

440  X  14 

0 

Tun- to 

0414 

141143 

37.24 

116.43 

ft 

65 

31 

5.4 

N 

Ooryeu 

0425 

1*3800 

970 

36.89 

11S.94 

T 

4 

4.5 

Pin  5lrip« 

0504 

133217 

646 

37.14 

116.14 

-A 

1.09  204 

190  X  17 

0 

Traveler  j 

osos 

140000 

1001 

37.05 

Uia.04 

A 

1.6  205 

548  X  55 

13 

8 

4.4 

ft 

Cycljnen 

0512 

193726 

aio 

37.13 

116.07 

A 

4.14  E04 

300  X  24 

10 

4.3 

A 

Tapestry 

0513 

133000 

1800 

37.09 

116.03 

T 

1.1  206 

1196  X  83 

100 

S.6 

ft 

Pir*t,hj 

0519 

135628 

2200 

37.11 

116.06 

T 

2.03  206 

1200  X  105 

190 

5.9 

N 

Diatom 

0527 

200000 

1106 

37.18 

116.10 

T 

3.9  205 

954  X  67 

21 

17 

5.0 

ft 

Discus 

Thrower 

0802 

153000 

1518 

37.23 

116.06 

Granite 

5* 

5.6 

Pile 

Ori*er 

0603 

140000 

1839 

37.07 

116.03 

T 

1.1  206 

1362  X  69 

140 

5.7 

N 

Tan  ! 

0615 

180247 

1494 

37.17 

116.05 

0 

7.01  205 

1300  X  60 

0 

kwiikuke* 

0625 

171300 

1057 

37.15 

116.07 

A 

2.41  205 

526  X  77 

25 

ft 

VvKjii 

0630 

221500 

2688 

37.32 

116.30 

ft 

1300  X  35 

300 

450 

6.1 

ft 

Hal  fU?4k 

0912 

ISjOOI 

835 

26.28 

115.95 

A 

12 

4.6 

0 

Cerriihjer 

0929 

144530 

750 

37.17 

116.05 

A 

8.54  204 

264  X  10 

4 

N 

?«la*<7<  1 

1105 

144500 

650 

37.17 

116.35 

A 

1.36  204 

190  X  15 

0 

Sit.-i 

1111 

120000 

782 

37.13 

116.05 

A 

6.33  204 

400  X  45 

0 

Ajti* 

ina 

1S0200 

693 

37.04 

116.01 

A 

9.9  204 

452  X  50 

0 

Cei  iso 

1213 

210000 

800 

36.88 

115.94 

A 

4.74  204 

200  X  125 

10 

4.6 

A 

*»i*w  Pc  mt 

1220 

153000 

825 

37.30 

116.41 

T 

(cyl Indrl- 

170  X  135 

825 

830 

6.3 

A 

Gnu  Icy 

cal ) 

1967 

0119 

n«500 

1194 

37.14 

116.13 

Ooloal te 

8.52  205 

500  X  180 

49 

5.3 

rtjin 

0120 

174003 

1836 

37.16 

116.00 

LI  rattan* 

35  X  135 

29 

5.3 

UOu'ltan 

0203 

151500 

844 

37.17 

116. 05 

A 

9.07  204 

260  X  30 

10 

4.6 

A 

Ward 

0223 

188400 

981 

37.02 

116.02 

A 

3.7  204 

560  X  20 

3 

4.4 

A 

Peis  i«m» 

0223 

185000 

2400 

37.13 

116.07 

T 

7.12  £05 

900  X  40 

130 

5.6 

ft 

A-Ji \e 

0302 

150000 

890 

37.17 

116.05 

A 

7.24  £04 

460  X  12 

0 

8l.tr  III 

0407 

150000 

889 

37.05 

116.02 

A 

1.4  205 

510  X  52 

0 

f  *wft 

0421 

150900 

789 

37.02 

116.06 

A 

4.2  204 

400  X  33 

7 

A 

Chocolate 

0427 

144500 

719 

37.14 

116.06 

A 

1.34  £04 

114  X  12 

0 

tr fundi 

0510 

134000 

1639 

37.08 

115.99 

T 

9  £33 

184  I  22 

10 

4.9 

A- 

NlCkCy 

0520 

150000 

2449 

37.13 

116.06 

T 

1120  X  148 

250 

230 

5.8 

Coerce  one 

0523 

140000 

3207 

37.27 

116.37 

T 

150 

NO 

5.7 

Stolen 

C526 

150002 

2069 

37.25 

116.48 

ft 

n 

47 

5.4 

Knuker* 

bOCkrr 

1 


I 


Table  VI 
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Y  Ml* 

Ott* 

Shat  TIM 

Device 

Oe pu  (ft) 

Epicenter 

Latitude  (°N)  Lsnotltud*  (®W) 

Htdlui  « 

1971 

1008 

143000 

1240 

37.11 

116.04 

T 

COflt. 

1106 

220000 

5875 

51.47 

-179.11 

B*  wit 

1214 

210959 

1085 

37.12 

116.09 

A 

1972 

0517 

14100 

1059 

37.12 

116.09 

A 

0720 

171600 

1391 

37.21 

116.18 

T 

0921 

153000 

1838 

37.08 

116.04 

T 

0926 

143000 

970 

37. 12 

116.09 

A 

1221 

201500 

2258 

37.14 

116.08 

T 

1973 

0308 

161000 

1866 

37.10 

116.03 

T 

0425 

222500 

1486 

37.00 

116.03 

A 

0426 

171500 

1850 

37.12 

116.06 

T 

0517 

160000 

39.79 

108.37 

Stndstone 

0605 

170000 

1284 

37.18 

116.21 

T 

0606 

130000 

3490 

37.24 

116.35 

T 

0628 

191512 

1530 

37.15 

116.09 

A 

1012 

170000 

1350 

37.20 

116.20 

T 

1974 

0227 

170000 

37.10 

116.05 

0619 

160000 

37.20 

116.19 

0710 

160000 

37.07 

116.03 

0830 

150000 

37.15 

116.08 

0926 

150500 

37. 13 

116.07 

1975 

0228 

151500 

37.11 

116.06 

0307 

150000 

37.13 

116.08 

0405 

194500 

37.19 

116.21 

0424 

141000 

37.12 

116.09 

0514 

14 0000 

2510 

37.22 

116.47 

0603 

142000 

2398 

37.34 

116.52 

0603 

144000 

2090 

37.09 

116.04 

0619 

130000 

2992 

37.35 

116.32 

0626 

123000 

4301 

37.28 

116.37 

1024 

171126 

440 

37.22 

116.18 

1028 

143000 

4150 

.37.29 

116.41 

1120 

150000 

2680 

37.22 

116.37 

1220 

200000 

2349 

37.13 

116.06 

1976 

0103 

191500 

4761 

37.30 

116.33 

0204 

142000 

2100 

37.07 

116.03 

0204 

144000 

2149 

37.11 

116.04 

0212 

144500 

3999 

37.21 

116.49 

0214 

113000 

3229 

37.24 

116.42 

0309 

140000 

2851 

37.31 

116.36 

0314 

123000 

4177 

37.31 

116.47 

0317 

141500 

2884 

37.26 

116.31 

0317 

144500 

2559 

37.11 

116.05 

Olu  I  Ht.  (ft) 

Ameuficdd 

CsttMUd 

\  use) 

» tH 

840  I  26 

7 

Cathey 

3600  I  55 

5000 

6.6 

CMAikiR 

372  X  41 

24 

Cheenjctl: 

336  X  48 

8 

Zimie 

21 

4.9 

Diamond 

>cu 1 1  & 

1350  X  90 

130 

5.6 

Oscuro 

388  X  51 

15 

15 

4.1 

Delpnioiu. 

584  X  97 

27 

4.8 

File 

1116  X  41 

67 

5.3 

Mi  ere 

21 

4.5 

Angus 

1150  X  125 

85 

120 

5.6 

Stjrwort 

90 

5.1 

*io  flllAC 

26 

5.0 

Otdu  Qute 

570 

6.1 

Alictnlni 

60 

4.9 

Portuljci  . 

9 

4.7 

Hutky  Act 

ISO 

5.6 

Ulfr 

20 

4.8 

King  Bled 

170 

5.7 

Cscabosg 

200 

5.6 

Poruun- 

ttfJU 

100 

5.5 

S tiny on 

IBS 

5.6 

TopijiUen 

120 

5.4 

Cibrillo 

20 

4.9 

Dining  Ce 

9 

4.5 

Etluh 

3  BO 

5.0 

T  ybo 

275 

5.B 

Stilton  : 

160 

5.6 

Mi22*n 

520 

5.9 

Musi 

750 

6.1 

15 

4.7 

Husky  P14 

1200 

8.2 

500 

5.9 

Inlet 

160 

5.6 

Cmberu 

600 

6.2 

Punster 

200 

5.6 

reel  son 

150 

5.6 

Cirwfl 

900 

6.1 

F online 

350 

5.8 

Chesnire 

350 

5.8 

Cstuui'jr 

900 

6.2 

Colby 

500 

6.0 

POvl 

200 

5.8 

Strait 

'  * 
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Origin 

H»r 

0«U 

Tin# 

Latitude  (°K) 

Longitude  (°E) 

(Bolt's) 

Announced 

tsttneted 

n^  (ISC)  Lscetton 

1964 

0315 

75958 

49.70 

78.00 

6.2 

49 

5.6 

0516 

60058 

49.90 

78.30 

6.2 

44 

5.6 

0719 

55959 

49.90 

78.10 

6.0 

29 

5.4 

0918 

75955 

72.90 

55.20 

2 

4.2  N.Z* 

1025 

75959 

73.50 

53.70 

5.3 

14 

5.1  N.Z. 

1116 

55957 

49.70 

78.00 

6.1 

49 

5.6 

ms 

0I1S 

55959 

49.89 

78.97 

7.0 

125 

no 

5.8 

0303 

61457 

49.82 

78.07 

6.0 

34 

S.S 

0S11 

63958 

49.79 

77.92 

6 

4.9 

0617 

24458 

49.97 

78.07 

5.8 

21 

S.Z 

0917 

3S9S8 

49.81 

78.05 

5.5 

15 

5.2 

1008 

55959 

49.89 

78.05 

5.8 

34 

5.4 

1121 

45758 

49.77 

78.06 

6.1 

47 

5.6 

1224 

45958 

49.88 

78.04 

7 

5.0 

1966 

0213 

45758 

49.82 

78.13 

6.5 

270 

6.1 

0320 

54958 

49.70 

78.00 

170 

6.0 

0421 

35758 

49.81 

78.05 

5.3 

28 

5.3 

0507 

35758 

49.74 

77.90 

4 

4.8 

0629 

65758 

49.93 

78.01 

36 

5.6 

0721 

35758 

49.70 

78.00 

5.9 

24 

5.3 

0805 

35758 

49.90 

78.00 

6.1 

29 

5.4 

0819 

35301 

50.40 

77.90 

4.6 

4 

S.l 

0907 

35158 

49.90 

78.00 

4 

4.8 

0930 

65953 

38.80 

64.50 

5.3 

30 

Uzbekistan  ■ 

1019 

35758 

49.75 

78.03 

6.3 

65 

5.6 

1027 

SS758 

73.44 

54.75 

6.5 

770 

8.4  N.Z. 

1218 

45758 

49.93 

77.73 

6.5 

120 

S.8 

1967 

0226 

35758 

49.78 

78.12 

6.6 

210 

8.0 

0325 

55759 

49.77 

78.08 

5.9 

21 

5.3 

0420 

40758 

49.74 

78.12 

6.3 

58 

5.5 

0528 

40758 

49.81 

78.11 

6.2 

32 

5.4 

0629 

25658 

49.87 

78.10 

27 

5.3 

0715 

32657 

49.83 

78.11 

6.Q 

27 

5.4 

0804 

65758 

49.82 

78.05 

5.8 

23 

5.3 

0916 

40358 

50.01 

77.82 

6.0 

22 

5.3 

0922 

50358 

50.03 

77.61 

6.0 

16 

5.2 

1017 

50358 

49.82 

78.10 

6.1 

67 

5.6 

1021 

45958 

73.37 

54.81 

6.0 

210 

5.9  N.2. 

1030 

60358 

49.84 

78.11 

6.0 

33 

5.3 

1(22 

40357 

49.90 

77.3 0 

3 

4.8 

1208 

60357 

49.84 

78.22 

22 

5.4 

1968 

0107 

34658 

49.81 

78.02 

10 

5.1 

0424 

103557 

49.83 

78.08 

7 

5.0 

0611 

30558 

49.84 

78.16 

5.8 

18 

5.2 

0619 

50557 

49.96 

79.09 

6.5 

35 

5.4 

0701 

40202 

47.92 

47.95 

5.7 

46 

5.S  N.  Caspfan  Sc# 

0712 

120757 

49.67 

78.12 

5.9 

23 

5.3 

082 Q 

40558 

50.00 

78.00 

4 

4.8 

0905 

40557 

49.7$ 

78.14 

6.2 

35 

5.4 

0929 

34258 

49.77 

78.19 

6.3 

110 

5.8 

1107 

100205 

73.40 

54.86 

6.0 

310 

8.1  N.2. 

1109 

25358 

49.79 

78.04 

4 

4.9 

1218 

50157 

49.72 

78.06 

5.7 

14 

1969 

0307 

82658 

49.81 

78.15 

6.3 

47 

5.6 

0516 

40257 

49.77 

78.15 

6.0 

18 

5.2 

0531 

50157 

49.98 

77.73 

6.2 

25 

S.3 

0704 

24657 

49.75 

78.19 

6.0 

22 

5.2 

0723 

24658 

49. 87 

78.32 

6.1 

38 

5.4 

0902 

45957 

57.41 

54.86 

5.2 

8 

11 

4.8  Urals 

0908 

45956 

57.36 

55.11 

5.2 

8 

11 

4.6  Ur*  Is 

0926 

65956 

4S.89 

42.47 

5.4 

78 

5.6  N.  Caspian  Sc* 

1001 

40258 

49.81 

78.21 

5.9 

21 

5.2 

1014 

70006 

73.40 

54.31 

6.5 

340 

6.3  N.Z. 

1130 

13257 

49.92 

79.00 

6.9 

160 

6.0 

1206 

70257 

43.83 

54.78 

5.7 

100 

5.8  6.  Caspian  V* 

1228 

34658 

50.00 

77.82 

6.5 

72 

5.7 

1229 

40158 

49.73 

78.15 

\ 

5.1 

1970 

0129 

70258 

49.80 

78.21 

S.9 

52 

5.5 

0327 

50257 

49.76 

78.01 

5.4 

10 

5.0 

0625 

4S9S2 

52.20 

55.69 

5.3 

6 

Ur*lt 

0628 

15758 

49.83 

78.25 

6.2 

120 

5.7 

0721 

30257 

49.95 

77.75 

6.0 

29 

5.4 

0724 

35657 

49.80 

78.17 

5.8 

21 

5.3 

0906 

40257 

49.77 

78.09 

6.0 

46 

5.4 

1014 

55957 

73.31 

55.15 

6.7 

6000 

2100 

5.6  N.Z. 

1104 

60257 

49.97 

77.79 

6.0 

34 

5.4 

1212 

70057 

43.85 

54.77 

6.6 

190 

6.0 

1217 

70057 

49.73 

78.13 

6.1 

35 

S.4 

1223 

70057 

43.83 

54.85 

6.6 

240 

6.0  £.  Caspian  S»* 

1971 

0322 

43258 

49.74 

78.18 

6.0 

26 

5.7 

0323 

65956 

61.29 

56.47 

5.9 

45 

51 

5.5  Urals 

0425 

332  58 

49.82 

78.09 

6.4 

140 

5.9 

0606 

40257 

49.98 

77.77 

5.5 

39 

5.5 

0619 

40358 

50.01 

77.74 

5.4 

36 

5.4 
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osu 

35057 

49.97 

79.05 

5.9 

0710 

165959 

04.17 

55.18 

5.1 

0919 

110007 

57.78 

41.10 

0927 

55955 

73.39 

55.10 

1009 

00257 

50.00 

77.70 

1021 

00257 

49.99 

77.65 

1022 

50000 

51.57 

54.54 

1129 

00257 

49.70 

78.13 

1215 

75259 

49.98 

77.90 

1222 

05950 

47.87 

48.22 

1230 

02058 

49.75 

78.13 

0210 

50257 

49.99 

78.89 

0.3 

0310 

45057 

49.75 

78.18  ' 

5.8 

0320 

42157 

49.73 

78.19 

5.6 

0411 

60005 

37.37 

62.00 

4.8 

0007 

I27S7 

49.70 

78.17 

5.7 

0700 

10258 

•49.72 

77.98 

4.8 

0709 

05958 

49.78 

35.40 

5.0 

0714 

145949 

50.00 

46.40 

3.6 

0310 

31057 

49.76 

78.15 

5.6 

0820 

25958 

49.40 

48.18 

6.3 

0820 

34057 

49.99 

77.78 

5.8 

0320 

55957 

73.34 

55.08 

0902 

85058 

49.96 

77.73 

5.3 

0921 

90001 

52.13 

51.99 

5.2 

1003 

85958 

46.85 

45.01 

6.1 

1102 

12058 

49.91 

78.84 

1124 

90008 

52.78 

51.07 

5.1 

1124 

95958 

51.84 

04.15 

5.1 

1210 

42058 

49.85 

78.18 

6.0 

1210 

42708 

50.11 

78.81 

6.7 

1220 

42713 

51.70 

79.20 

4.9 

0210 

50258 

49.83 

78.23 

5.5 

0419 

43258 

50.01 

79.72 

0710 

12658 

49.78 

78.06 

0723 

12258 

49.99 

78.85 

0315 

15958 

42.71 

67.41 

0828 

25958 

50.55 

08.39 

0912 

65954 

73.30 

55.16 

0919 

25957 

45.63 

07.85 

0927 

65958 

70.76 

53.87 

0930 

45957 

51.61 

54.58 

1020 

42658 

49.70 

78.20 

1020 

55958 

53.66 

55.38 

1027 

65957 

70.78 

54.18 

1214 

74657 

50.04 

79.01 

0130 

45658 

49.89 

77.99 

0130 

45702 

49.83 

78.08 

0410 

55302 

49.99 

78.82 

0516 

3025 7 

49.74 

78.15 

0S31 

32657 

49.95 

78.84 

0625 

35658 

49.89 

78.11 

0710 

25057 

49.79 

78.14 

0814 

1459a 

68.91 

75.90 

0829 

959» 

73.37 

55.09 

0829 

150000 

67.23 

62.12 

0913 

30258 

49.82 

76.09 

1016 

63257 

49.97 

78.97 

1102 

45957 

70.  S2 

54.06 

1207 

55957 

49.91 

77.65 

1216 

62302 

49.75 

78.06 

1210 

64102 

49.82 

78.12 

1227 

54657 

49. 96 

79.05 

0220 

53258 

49.82 

78.08 

0311 

542a 

49.79 

78.25 

0427 

53057 

49.99 

78.98 

0008 

320a 

49.76 

78.09 

0807 

35658 

49.81 

78.24 

0823 

85958 

73.37 

54.04 

0929 

105958 

69.51 

90.40 

1018 

85956 

70.84 

53.69 

1021 

115957 

73.35 

55.08 

1029 

44658 

49.98 

78.97 

1213 

4US7 

49.80 

78.20 

1225 

51657  . 

50.04 

78.90 

0115 

44658 

49.87 

78.25 

0421 

50257 

49.93 

78.82 

0704 

zasi 

49.91 

78.95 

0723 

232a 

48.79 

78.05 

25 

5.2 

27 

5.2 

Uralt 

4 

4.5 

Urili 

770 

0.5 

8.2. 

24 

5.3 

44 

5.5 

34 

5.2 

Uralt 

34 

5.4 

3 

4.9 

210 

0.0 

8.  Caspian  5«a 

90 

5.7 

43 

5.4 

33 

5.4 

15 

5.1 

7 

Turtaaa 

34 

5.4 

1 

4.4 

6 

4.8 

N.  Slack  Sea 

0.2 

8.  Caspian  Sm 

15 

5.0 

87 

5.7 

8.  Caspian  Sna 

3S 

5.3 

090 

0.3 

8.2. 

7 

4.9 

21 

8.  Caspian  5m 

88 

5.0 

IM  Caspian  5c a 

350 

0.1 

11 

4.5 

Urals 

20 

5.2 

W.  Kazakh 

70 

5.0 

020 

0.0 

28 

5.3 

Uzbekistan 

14 

5.2 

u. 

Kazakh 

2700 

6.8- 

8.2. 

11 

5.1 

u. 

Kaiakh 

210 

5.9 

8.2. 

22 

5.2 

Urals 

19 

5.2 

7 

4.8 

Uralt 

3200 

6.9 

8.2. 

150 

5.8 

2 

4.9 

23 

5.2 

3 

4.9 

23  - 

5.2 

140 

5.9 

2 

4.7 

10 

5.2 

45 

5.4 

8.2. 

870 

6.4 

8.2. 

20 

5.0 

Uralt 

15 

5.2 

43 

5.5 

1600 

0.4 

8.2. 

2 

4.1 

a 

5.0 

0 

4.8 

51 

5.6 

77 

5.7 

30 

5.4 

60 

5.6 

35 

5.5 

14 

5.2 

550 

6.3 

8.2. 

6 

4.8 

u. 

.  Siberia 

1400 

6.7 

8.2. 

700  ■ 

6.6 

8.2. 

90 

5.8 

10 

S.l 

90 

5.7 

14 

5.2 

20 

5.3 

90 

5.8 

10 

5.1 

*  r 
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